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本 书 从 实际 应 用 的 角度 ， 研 究 并 改进 了 现 有 的 图 像 检测 与 分 割 方法 。 
重点 讲述 了 基于 圆 点 滤波 器 的 增强 算法 、 组 合 优化 〈 如 改进 遗传 算法 或 分 
BR) 和 支持 向 量 机 相 结 合 的 “去 假 阳 ”算法 框架 、 基 于 主动 形状 模 
型 的 分 割 方法 、 基 于 多 尺度 的 均值 漂移 的 带宽 选取 方法 、 结 合 均值 漂移 和 
期 望 最 大 的 分 割 方法 等 。 

本 书 以 肺 部 的 孤立 型 肺 结 节 图 像 检 测 和 分 割 为 实际 案例 ， 在 应 用 中 对 
这 些 算 法 进行 了 改进 研究 ， 且 提出 了 结合 均值 漂移 和 期 望 最 大 的 图 像 分 割 
方法 ; 通过 在 肺 部 图 像 的 粘连 血管 型 肺 结 节 分 割 中 的 应 用 ， 证 明了 其 具有 
较 好 的 应 用 效果 。 还 提出 改进 主动 形状 模型 方法 ， 并 应 用 在 肺 部 图 像 的 病 
态 肺 实质 分 割 中 ， 取 得 了 和 较 好 的 结果 。 

本 书 详细 介绍 了 这 些 算法 的 理论 推导 、 算 法 步骤 和 实验 结果 ,力求 向 
读者 展示 出 图 像 检 测 与 分 割 的 最 新 研究 动态 及 其 应 用 成 果 ， 和 希望 能 为 从 事 
图 像 处 理 的 广大 读者 提供 参考 ， 推 动 该 领域 的 学 术 和 进步。 本 书 理论 联系 实 
际 ， 应 用 实例 清晰 可 靠 。 

本 书 第 1、2 章 由 任 会 之 编写 ,第 3、4、5 章 由 孙 申 申 编写 ， 全 书 由 
任 会 之 统 稿 。 

本 书 的 主要 内 容 为 著者 多 年 来 的 研究 成 果 。 本 书 的 研究 课题 得 到 了 国 
家 自然 科学 基金 青年 基金 (71201105) 、 辽 宁 省 高 校 科学 研究 项 目 
(12014480) 的 资助 。 

由 于 图 像 检测 与 分 割 技术 正在 迅速 发 展 之 中 ， 加 之 作者 水 平 有 限 ， 本 
书 的 内 容 会 有 不 足 之 处 ， 奶 请 广大 读者 批评 指正 。 
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$1% JET Hessian 年 阵 的 球形 
区 域 检 测 方 法 及 其 应 用 


所 谓 图 像 检测 ， 就 是 对 图 像 感 兴趣 的 特征 区 域 (检测 目标 ) 有 针对 性 
地 进行 提取 及 其 过 程 。 其 中 图 像 是 承载 检测 目标 的 载体 ， 检 测 目标 需要 事 
先进 行 特征 提取 、 上 归纳， 最 终 通 过 相应 算法 分 离 出 来 。 本 章 试 图 应 用 基于 
Hessian 矩阵 的 圆 点 滤波 融 ， 对 具有 球形 特征 的 目标 图 像 进行 检测 。 

















1.1 基于 二 维 Hessian 和 矩阵 的 圆 点 滤波 器 


























在 二 维 上 ， 服 从 高 斯 分 布 的 圆 形 用 式 (1-1) 表示 ， 服 从 高 斯 分 布 的 
线 结构 用 式 (12) RR, 





dla) =op =) GD) 
(x,y) =exp[ 2] (12) 

二 维 Hessian 矩阵 形式 如 下 : 
=) 本 


村 征 值 为 A,、A ， 当 和 =A, <0 时 ， 为 圆 上 的 像素 ; 当 和 ,<0 HA, = 
0 时 ， 为 线 上 的 像素 。 在 二 维 空间 下 ， 定 义 的 圆 点 和 线 的 增强 算 子 如 下 : 
[A,?/|A,|, A, <0,A, <0,A, <0 时 


À, À, À) = 1-4 
Ban ( 了 ii 3) y 其 他 ( ) 


IA,1-1A,1, “4A, <0 Mt 
| 1 2 1 (1-5) 


Ziine (A1 A2,A3) = 0 其 他 


1.2 基于 三 维 Hessian 矩阵 的 圆 点 滤波 器 











在 三 维 上 ， 用 式 (1-6) 表示 服从 高 斯 分 布 的 球 ; 用 式 (1-7) 表示 服 
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从 高 斯 分 布 圆 柱 ; 用 式 (1-8) 表示 服从 高 斯 分 布 的 面 模型 。 模 型 方程 
wW Fi”: 


— (x? a +2) 





daina = exp| = (1-6) 
Kuna) sap A] (1-7) 
r (8) 
三 维 Hessian 和 矩阵 形式 如 下 : 
APE (1) 





REED A, AD. Ag, 4A, =A, =A, <0 时 ， 为 结 节 上 的 像素 ; 当 
A, =A, <0 H à, =0 时 ， 为 血管 上 的 像素 ; 当 A, <0 Hà, =A, =0 时 ， 为 
平面 上 的 像素 。 

图 1-1 阐述 了 Hessian 和 矩阵 特征 值 和 局 部 形状 之 间 的 关系 。 在 三 维 空 
间 中 ， 和 定义 的 管 、 平 面 和 球 的 增强 算 子 如 下 : 





a) b) c) 
图 1-1 Hessian 和 矩阵 特征 值 和 局 部 形状 之 间 的 关系 
a) 管 b) 平面 c) 球 











à PIA |, HA, <0,A, <0,A, <0 时 
KEETE sl/l ita ° . (1-10) 
A à l- la, DZIà l, “4A, <0,A, <0,A, <0 时 

zm Ar Asda) = 2 [C | | ;1) | | ith 7 3 
(1-11) 


(1-12) 


IA l-11, 4A, <0 FY 
Zilane (A1 »A2 ,A ) = 


0, 其 他 
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1.3 圆 点 滤波 器 理论 在 图 像 处 理 中 的 应 用 实例 


圆 点 滤波 需 理论 可 以 应 用 在 具有 球形 特征 的 图 像 目标 检测 中 。 在 特定 
人 群 的 人 体 肺 部 CT 影像 上 ， 毛 玻璃 型 肺 结 节 (CGO) 就 表现 为 球形 特 
征 。 而 检测 这 些 存在 于 CT 影像 上 的 毛 玻璃 型 肺 结 节 ， 有 一 定 的 临床 价值 。 
因为 圆 点 滤波 器 可 以 增强 球形 图 像 特征 5 ， 所 以 用 圆 点 滤波 器 增强 毛 玻 璃 
肺 结 节 并 进行 检测 是 可 行 的 。 其 具体 方法 和 步骤 如 下 。 

首先 ， 利用 自 适 应 非 线 性 滤波 器 (AN filter, Adaptive Nonlinear fil- 
ter) MARE SERN AA LEE; PR, A e PA 
把 血管 从 肺 实质 区 中 噜 除 掉 ， 这 时 ， 肺 实质 区 内 只 有 毛 玻 璃 型 结 节 和 肺 实 
质 ; 最 后 ， 用 圆 点 滤波 器 就 可 以 提取 出 候选 毛 玻 璃 型 结 节 。 

之 所 以 用 自 适应 非 线性 滤波 需 是 因为 它 是 空域 上 的 滤波 器 ， 不 用 做 传 
里 叶 变 换 :“ ， 速 度 上 是 可 以 接受 的 。 同 时 ， 由 于 其 自 适应 性 ， 它 在 去 除 噪 
声 的 同时 ， 也 增强 了 边缘 ， 且 有 较 好 的 拉 伸 对 比 度 效 果 。 检 测 GGO 结 节 
的 算法 程序 流程 图 如 图 1-2 所 示 。 


输入 影像 数据 
自 适 应 滤波 器 
二 值 化 去 除 血管 


二 值 化 提取 高 亮点 



































圆 点 滤波 器 提取 
GGO “候选 点 ” 





图 12 GGO 结 节 检 测算 法 程序 流程 图 








由 于 固定 模板 的 缺点 ，Guillon 等 人 ”提出 了 根据 图 像 的 局 部 梯度 自 
适应 地 变换 滤波 器 模板 ， 即 自 适 应 非 线性 滤波 器 (AN filter) 。 设 模板 S 的 
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尺寸 为 K,. K, ， 模 板 的 系数 应 该 是 像素 (x,y) 和 像素 (x +i,y +j) RE 
值 相似 程度 。 如 果 相 同 ， 则 my = 1; 如 果 相 差 很 大 ， 则 my =0。 自 适应 滤 
波 器 的 模板 M= {m;e [0,1] (i,j) e S| ， 定 义 自 适应 滤波 器 模板 的 系数 
ms 如 下 : 





el eas =), Ez ia =i Aaa 


[f(x tisy +7) -f(x,y)] | 
28° 





m; =exp| - (1-14) 


H T EE Ar ERA RAER, SER a ES at 
和 高 通 滤波 器 结合 起 来 。 


g(x,y) =g,(%,y) +0g,(«,y) (1-15) 
g(x,y) = Dh fle +i,y +j) (1-16) 
(i,j) es 
h; = it (1-17) 
My | 
(k,l eS 
en(xy) = > w fla ti,y +j) (1-18) 
2 
W; j = m;;- m; m m. . (1-19) 


= a t 

式 中 : gj(x,y) 为 低 通 滤波 器 ;gj (x,y) 为 高 通 滤 波 器 ;Card (S) 表示 集合 
S 中 元 素 的 个 数 。 

自 适应 滤波 器 低 通 部 分 和 高 通 部 分 的 效果 如 图 1-3 所 示 。 图 1-3a 是 对 

含 GGO 的 CT 影像 进行 自 适应 非 线性 滤波 的 低 通 部 分 的 效果 ， 网 1-3b 是 

对 含 GGO 的 CT 影像 进行 自 适 应 非 线性 滤波 的 高 通 部 分 的 效果 。 图 1-4 为 





SS 


a) b) 


图 1-3 自 适 应 滤波 器 低 通 部 分 和 高 通 部 分 的 效果 
a) 低 通 部 分 效果 b) 高 通 部 分 效果 
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1-4 自 适应 非 线性 滤波 器 滤波 的 增强 效果 





自 适应 非 线性 滤波 器 滤波 的 增强 效果 ， 可 以 看 出 ， 滤 波 器 拉 伸 了 血管 和 毛 
玻璃 型 结 节 的 对 比 度 ， 也 拉 近 了 毛 玻 璃 型 结 节 和 肺 实质 区 域 的 CT 对 比 度 。 
如 果 数 据 集合 服从 正 态 分 布 W(w,rc) ， 根 据 正 态 分 布 的 性 质 ，4c 能 
盖 正 态 分 布 曲 线 下 面积 的 0.95 倍 。 图 1-4 所 示 的 直方 图 如 图 1-5a 所 示 。 
假设 该 直方 图 服从 半 个 正 态 分 布 ， 血 管 等 的 高 亮 区 域 应 该 占 曲 线 下 面积 也 

(1 -0.95)/2 =0. 025 倍 。 





1000 











a) b) 


图 1-5 从 增强 图 像 中 剔除 血管 
a) 图 1-4 的 直方 图 b) 去 除 血管 后 的 图 像 











图 1-4 所 示 中 可 以 看 到 已 经 把 血管 和 毛 玻 璃 型 结 节 的 对 比 度 拉 大 了 ， 
下 面 用 自 适 应 阔 值 的 方法 再 把 血管 剔除 出 去 。 设 肺 实质 区 域 直方 图 的 最 大 
灰 度 级 数 为 maxValue， 最 小 灰 度 级 数 为 minValue， 肺 实质 区 域 的 像素 个 数 
为 n， 从 maxValue 到 minValue 累加 灰 度 级 对 应 的 像素 个 数 和 sum Threshold , 
当 sumThreshold > n x0.025 则 停止 ， 并 记录 当前 的 灰 度 值 ThresholdValue。 
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HEME r BBC ThresholdValue, FAA {EL r 把 血管 从 原 图 像 中 剔除 掉 ， 


如 图 


1-5b 所 示 。 





如 果 直 接 对 图 1-5b 做 圆 点 滤波 器 的 处 理 ， 运 算 量 会 非常 大 ， 且 耗 时 
多 ， 不 实用 。 则 需要 对 图 1-5b 取 高 亮 区 域 ， 在 高 亮 区域 中 进行 圆 点 滤波 
器 的 操作 。 用 自 适 应 阔 值 的 方法 提取 高 亮 区 域 ， 即 设 图 1-5a 直方 图 的 最 
大 灰 度 级 数 为 maxValuel ， 最 小 灰 度 级 数 为 minValuel ， 肺 实质 区 域 的 像素 
ŽU ni, AÀ maxValuel 到 minValuel 累加 灰 度 级 对 应 的 像素 个 数 和 
sumThresholdl ， 当 sumThreshold1 >n, x0.1 则 停止 ， 并 记录 当前 的 灰 度 级 
数 ThreshodTuluel。 自 适应 国 值 fi 被 设 为 ThresholdValuel 。 二 值 化 后 高 亮 
区 域 如 图 1-6a 所 示 ; 对 图 1-6a 进行 圆 点 滤波 顺 操 作 结 果 如 图 1-6b ras, 
图 中 圈定 的 为 候选 毛 玻 璃 型 结 节 。 





























1.4 


1. 4. 


影像 
一 附 


a) 





图 1-6 提取 毛 玻 璃 型 候选 结 节 的 结果 


D) 








a) 高 亮 区 域 b) 圆 点 滤波 器 提取 得 GCO 型 “候选 点 ” 


实验 结果 及 分 析 


1 实验 数据 





表 1-1 描述 了 本 章 中 使 用 的 影像 数据 来 源 和 数据 属性 。 数 据 来 源 于 
LIDC 数据 库 和 广州 医学 院 第 一 附属 医院 ， 每 套数 据 的 结 节 都 有 医生 手工 
标识 的 “ 金 标准 "， 即 结 节 的 位 置 。LIDC 数据 库 中 的 数据 是 扫描 部 分 的 肺 


， 而 广州 医学 院 第 一 附 











属 医院 的 数据 是 扫描 全 肺 影像 。 广 州 医学 院 第 








属 医院 的 数据 库 中 的 影像 ， 是 通过 东芝 16 排 CT 产生 的 ， 球 管 的 电压 


为 120kV， 电 流 为 130mA。 数 据 被 重建 成 5312 x512 大 小 的 影像 。 
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表 1-1 本 章 测试 算法 性 能 的 数据 库 









































项 目 LIDC 数据 库 广州 医学 院 第 一 附属 医院 

CT 影像 套数 / 套 23 64 

像素 单位 体积 0. 6mm x0. 6mm x0. 6mm 0.6mm x0. 6mm x2mm ~ 3mm 

平均 层 数 / 层 30 ~50 100 ~200 

“ 金 标准 ”中 结 节 总 数 / 个 23 226 
孤立 型 结 节 总 数 / 个 11 82 
粘连 血管 型 结 节 总 数 / 个 6 54 
粘连 肺 壁 型 结 节 总 数 / 个 2 58 
GGO 结 节 数 /个 4 32 











1.4.2 实验 环境 


所 有 实验 都 是 在 主 频 AMD2GHz、 内 存 1.5GB 计算 机 上 完成 的 ; 操作 
系统 是 Windows XP， 编 程 语言 为 VC ++6.0。 算 法 代码 的 开发 都 是 在 东软 
医学 影像 处 理 与 分 析 平 台 下 完成 的 ， 该 平台 是 基于 NIL 库 实现 的 ， 并 负责 
把 3D 影像 数据 读 入 到 内 存 中 ， 且 能 把 处 理 后 的 结果 显示 到 界面 上 。 


1.4.3 实验 方案 


本 章 实验 共 设 计 2 套 实验 方案 ， 具 体 方 案 如 下 。 

1. 实验 一 

目的 是 分 析 提 取 毛 玻璃 型 候选 肺 结 节 的 必要 性 。 让 医生 检测 LIDC 数 
据 库 中 的 所 有 影像 ， 并 画 出 检测 结果 。 同 时 和 LIDC 发 布 的 “ 金 标准 ” 进 
行 比较 ， 并 分 别 统计 各 种 类 型 结 节 的 漏 检 情况 。 

2. 实验 二 

目的 是 分 析 所 提出 的 提取 毛 玻 璃 型 候选 肺 结 节 方 法 的 有 效 性 和 实用 
性 。 首 先 ， 找 出 “人 金 标准 ”含有 毛 玻 璃 型 结 节 的 影像 ; 然后 ， 分 别 对 每 套 
影像 分 割 肺 实质 区 域 ; 再 次 ， 按 照 1.3 中 所 述 方法 ， 得 到 Coco 结 节 的 
“候选 点 ”; 最 后 ， 统 计 平 均 运 行 时 间 、“ 候 选 点 ”个 数 / 套 、 漏 检 的 个 数 。 
与 此 同时 ， 在 相同 的 数据 集 上 对 文献 [6，8] 所 述 算法 进行 了 评估 ， 并 把 
它们 的 结果 与 所 提出 的 方法 的 结果 相 比 较 。 


1.4.4 实验 结果 


实验 一 的 实验 结果 是 : LIDC 数据 库 中 共 23 个 结 节 ， 医 生 总 共 漏 检 了 6 
个 ， 其 中 孤立 型 结 节 汤 检 2 个 ,粘连 血管 型 结 节 和 粘连 肺 壁 型 结 节 无 漏 检 ， 而 
毛 玻 璃 型 结 节 漏 检 了 4 个 〈 漏 检 率 100% )。 因 此 ,证 明了 毛 玻璃 型 结 节 是 医 
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生 最 易 漏 检 的 结 节 类 型 ， 从 而 说 明 提 取 候 选 毛 玻 璃 型 结 节 的 必要 性 。 
实验 二 中 把 1. 3 所 述 方法 中 的 自 适应 滤波 器 的 参数 a 设 为 3， 模 板 的 
高 斯 系数 delta 也 设 定 为 3。 检 测 结果 是 : 对 于 广州 医学 院 第 一 附属 医院 数 
据 库 中 的 32 个 毛 玻 璃 型 结 节 ， 漏 检 了 5 个 ; 而 对 于 LIDC 数据 库 中 的 4 个 
毛 玻璃 型 结 节 ， 没 有 漏 检 现象 发 生 。 表 1-2 分 别 列 出 了 3 种 方法 的 漏 检 
率 、 平 均 每 层 “ 候 选 点 ”个 数 ， 平 均 每 套 影像 提取 “候选 点 ”所 需 时 间 。 


#12 本章 方 法 与 现 有 的 2 种 方法 的 性 能 对 比 


























实验 方法 漏 检 率 候选 点 个 数 / 层 速度 / (min/ 套 ) 
文献 [6] 方法 42. 5% 32.0 3.5 
文献 [8] 方法 30. 4% 7.4 4.5 

本 章 方 法 15.6% 5.5 1.2 











图 1-7 为 4 幅 包 含 GGO ZAHIR, SPAN GGO HW, AR 
方法 提取 候选 毛 玻璃 型 结 节 的 实际 效果 如 图 1-8 所 示 ， 图 1-8 为 本 章 算 法 
提取 GGO 结 节 “候选 点 ”( 各 分 图 的 左 半 部 分 ) 的 结果 ， 且 用 箭头 标示 。 





a) b) 





à 











到 1-7 4 幅 包 含 GGO 结 节 的 影像 
a) 含 毛 玻璃 型 结 节 影 像 1 ob) 含 毛 玻璃 型 结 节 影像 2 
c) 含 毛 玻璃 型 结 节 影 像 3 d) 含 毛 玻璃 型 结 节 影 像 4 
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c) 
图 1-8 本 章 算 法 提取 GGO 结 节 “候选 点 ”的 结果 


a) 对 影像 1 提取 “候选 点 ” 
c) 对 影像 3 提取 “候选 点 ” 
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b) 对 影像 2 提取 “候选 点 ” 
d) 对 影像 4 提取 “候选 点 ” 
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第 2 革 基于 组 合 优化 和 SVM 的 
图 像 球形 区 域 检测 方法 





在 模式 识别 和 图 像 区 域 检测 中 ,具有 上 自我 学 习 能 力 的 支持 向量 机 
(SVM, Support Vector Machine) 是 一 种 常用 方法 。 应 用 支持 向 量 机 ， 先 
要 对 图 像 进 行 特 征 提 取 。 在 提取 的 众多 特征 中 选取 部 分 特征 进行 组 合 ， 以 
实现 模式 识别 或 图 像 区 域 检测 的 最 优 效果 。 选 择 特 征 组 合 的 问题 可 以 看 作 
是 一 个 组 合 优化 问题 。 本 章 将 遗传 算法 、 分 散 搜索 算法 等 求解 组 合 优化 问 
题 的 方法 与 SVM 相 结合 ， 并 应 用 于 图 像 区 域 检测 。 


2.1 求解 组 合 优化 问题 的 智能 计算 方法 


























对 于 图 像 特征 选取 的 组 合 优化 问题 ,采用 了 求解 组 合 优化 问题 的 两 种 
现代 优化 算法 ， 即 遗传 算法 和 分 散 搜索 算 法 。 


2.1.1 遗传 算法 


WERE" (GA, Genetic Algorithm) 模拟 了 达尔 文 的 “优胜 劣 汰 ” 
进化 过 程 ， 即 模拟 了 自然 选择 和 和 遗传 进化 中 发 生 的 繁殖 、 交 配 和 突变 现 
象 。 从 一 个 初始 种 群 出 发 ， 通 过 随机 选择 、 交 叉 和 变异 操作 ， 产 生 一 群 新 
的 更 适合 环境 的 个 体 ， 再 使 群体 进化 到 搜索 空间 中 越 来 越 好 的 区 域 ， 这 样 
一 代 代 不 断 进化 ; 最 后 ， 收 敛 到 一 个 最 适应 环境 的 个 体 上 ， 求 得 问题 的 近 
优 解 。 

在 遗传 算法 中 ,需要 用 到 适应 度 函 数 、 选 择 算 子 、 交 义 算 子 和 变异 血 
子 等 概念 。 适 应 度 是 为 实际 问题 预先 定义 的 数学 度量 ， 用 于 测试 每 个 个 体 
的 优 劣 程度 。 选 择 算 子 保证 父 种 群 和 子 种 群 个 体 的 连续 性 ， 也 使 得 子 种 群 
继承 父 种 群 适应 度 函 数 高 的 个 体 。 交 叉 算 子 和 变异 算 子 用 在 种 群 进化 中 ， 
交叉 算 子 使 得 下 一 代 种 群 中 的 个 体 不 但 继承 上 代 种 群 中 父 体 的 基因 ， 而 且 
还 要 继承 母体 的 基因 ; 变异 操作 对 应 着 生物 界 的 基因 突变 ， 增 加 种 群 中 个 
体 的 多 样 性 。 

每 个 用 二 进 制 串 编码 的 个 体 叫 作 染 色 体 。 遗 传 算法 的 具体 步 又 如 下 : 
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1) 建立 初始 种 群 ， 由 PP 个 个 体 组 成 ， 每 个 个 体 为 d 位 基因 码 ， 计 算 
每 个 个 体 的 适应 度 (Fitness) PRBCE 

2) 选择 算 子 (如 轮 盘 法 、 锦 标 赛 选择 等 方式 ) : 从 当前 种 群 中 选择 
繁衍 到 子 种 群 的 个 体 。 

3) 从 2) 选择 的 个 体 经 过 交叉 算 子 (如 两 点 交叉、 均匀 交叉 等 方式 ) 

和 变异 算 子 (如 “位 ”突变 ) ， 产 生 后 代 种 群 。 

4) 从 3) 中 选择 前 P 个 适应 度 函 数值 最 大 的 个 体 作为 当前 种 群 。 重 
复 步骤 2) ~4) 直到 终止 条 件 被 满足 。 

步骤 中 的 终止 条 件 可 以 是 事先 定义 的 最 多 迭代 次 数 和 求解 精度 等 。 


2.1.2 分 散 搜 索 算法 


与 遗传 算法 类 似 ， 分 散 搜索 (SS, Scatter Search)'" 算法 也 是 一 种 基 
于 种 群 的 进化 算法 。 不 同 的 是 ,分 散 搜索 算法 在 种 群 的 选择 中 ， 不 仅 考 虑 
个 体 解 的 质量 ， 还 要 考虑 整个 种 群 的 分 散 性 。 通 过 控制 参考 集中 解 的 质量 
和 分 散 性 ， 使 得 搜索 过 程 具有 良好 的 广泛 性 ， 不 容易 陷 和 人 局 部 最 优 。 算 法 
包括 5 个 部 分 : 初始 种 群 的 产生 策略 、 子 集 产 生 策略 、 解 的 组 合 策略 、 解 
的 改进 策略 和 参考 集 更 新 策略 。 

初始 种 群 的 产生 策略 必须 能 产生 具有 和 良好 质量 和 分 散 性 的 初始 种 群 ， 
这 样 从 中 精 选 出 的 参考 解 集 ， 才 能 广泛 地 分 布 在 整个 解 空间 ， 且 具有 良好 
的 多 样 性 。 子 集 产生 策略 构造 若干 个 参考 集 的 子 集 ， 而 解 的 组 合 策略 分 别 
组 合 这 些 子 集中 的 解 产生 新 解 。 产 生 的 这 些 新 解 经 改进 策略 进行 改进 后 ， 
用 于 更 新 参考 集 。 分 散 搜索 算法 步骤 如 下 : 

1) 以 初始 种 群 产生 策略 建立 初始 种 群 。 

2) 使 用 参考 集 更 新 策略 从 初始 种 群 构造 参考 集 ， 其 过 程 (参考 集 包 
含 新 解 ) 如 下 : 

D 用 子 集 产生 策略 产生 参考 集 的 子 集 集合 NewSubsets。 

(2) 对 NewSubsets 中 的 各 个 子 集 分 别 使 用 组 合 策略 产生 新 解 ， 并 以 改 
进 策略 进行 改进 。 

© 用 改进 后 的 这 些 新 解 更 新 参考 集 。 


2.2 基于 组 合 优化 的 特征 选择 模型 





































































































特征 选择 方法 有 两 种 模式 : 过 滤 (filter) 模式 和 封装 (wrapper) 模 


-12- 图 像 检测 与 分 割 方法 及 其 应 用 








BOS 。 过 滤 模 式 利 用 特征 对 数据 的 可 分 度 和 特征 间 的 相关 性 作为 选取 特 
征 组 合 的 指标 ， 一 般 具 有 效率 高 而 效果 差 的 特点 ; 封装 模式 利用 分 类 顺 的 
分 类 正确 率 作为 选取 特征 组 合 的 指标 ， 具 有 效率 低 但 效果 好 的 特点 。 由 于 
建立 分 类 器 是 在 离线 状态 下 进行 的 ， 对 速度 的 要 求 不 高 ， 而 且 该 分 类 器 用 
于 医学 背景 下 ， 需 要 保证 较 高 的 敏感 性 和 特异 性 ， 所 以 选择 封装 模式 较为 
适宜 。 

特征 选择 实际 是 一 个 特征 的 组 合 优化 问题 ， 已 被 证 明 是 NP-Hard 问 
题 ， 和 帝 见 的 智能 优化 方法 如 GA 被 应 用 到 解决 该 问题 中 。 对 于 许多 图 像 检 
测 的 应 用 (如 医学 影像 中 的 模式 识别 ) 上 的 标准 数据 ， 基 于 SVM 建立 的 
分 类 器 效果 要 优 于 人 工 神经 元 网 络 (ANN, Artificial Neural Networks) t°! 、 
决策 树 (DT, Decision Trees) 等 其 他 分 类 器 ,这 是 因为 SVM 是 专门 针 
对 有 限 样 本 分 类 情况 的 。SVM 最 终 把 分 类 问题 转化 成 为 一 个 二 次 型 寻 优 问 
题 ， 可 以 得 到 全 局 最 优 解 ， 解 决 了 人 工 神经 元 网 络 中 无 法 避免 的 局 部 极 值 
问题 。SVM 将 原 特征 空间 通过 非 线 性 变化 转换 到 高 维特 征 空间 ， 在 高 维特 
征 空间 中 构造 线性 判别 函数 来 实现 原 空 间 中 的 非 线 性 判别 函数 。 图 2-1 描 
述 了 基于 组 合 优 化 和 SVM 的 特征 选择 的 过 程 框架 。 


3 a 
训练 集 进行 特征 选择 
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候选 特征 组 合集 





初始 特征 组 合 


o> BS ES ae 









训练 SVM 改进 GA 
分 类 器 算法 评估 


图 2-1 基于 组 合 优化 和 SVM 的 特征 选择 的 过 程 框架 











在 图 2-1 中 首先 由 组 合 优 化 算法 产生 初始 候选 特征 组 合集 (初始 候选 
特征 子 集 ) ， 并 和 训练 集 数据 一 起 训练 出 SVM 分 类 器 ; 然后 ， 用 SVM 分 
类 器 测试 用 于 特征 选择 的 测试 集 ， 得 出 的 测试 结果 返回 给 组 合 优化 算法 进 
行 对 候选 特征 组 合 的 评估 ， 并 产生 下 一 代 候 选 特征 子 集 ， 如 此 反复 迭代 到 
最 大 代数 ， 从 而 求 出 最 优 特 征 组 合 。 对 特征 选择 问题 建立 数学 模型 如 下 : 

maxf( X) =w, » Sensitivity( X) +w, + Specificity(X ) (2-1) 
s.t.  Sensitivity(X) =TP(X)/[ TP(X) +TN(X)] 
Specificity(X) =FN(X)/[ FP(X) + FN(X) | 
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X,e {0,1} 
式 中 : w +w, =1 并 且 w >w,; Sensitivity 为 敏感 性 率 ; Specificity 为 特异 性 
率 ， 检 测 的 敏感 性 和 特异 性 之 和 为 1; TP 为 真 阳性 ; FP 为 假 阳 性 ; FN 为 
(BAKE; TN 为 真 阴 性 ; w 和 w, 分 别 为 敏感 性 和 特异 性 的 权重 参数 ， 它 们 
是 [0,1] 之 间 的 实数 ; /为 [0,1] 之 间 的 实数 ; X 为 特征 组 合 。 


2.3 基于 改进 GA 和 SVM 的 求解 模型 





首先 ， 定 义 染色 体 ;， 染色 体 为 0-1 二 进 制 串 ， 其 长 度 为 所 提取 特征 序 
ATR, 1 表示 该 特征 被 选中 ， 否 则 为 没有 被 选中 。 其 次 ， 设 定 初 始 
种 群 ， 随 机 生成 长 度 为 R 的 染色 体 w 个， 作为 初始 种 群 的 个 体 ， 为 了 保 
证 种 群 中 个 体 的 多 样 性 ， 本 例 中 取 N =80。 再 次 ,设计 适应 度 函 数 ; 适应 
ERR SVM 训练 出 的 分 类 融 的 性 能 有 关 ， 其 中 包括 敏感 性 和 特异 性 ; 
定义 适应 度 函 数 为 式 (2-1) 模型 中 的 目标 函数 。 在 分 类 效果 相同 的 情况 
下 ， 竺 征 组 合 中 所 含 特征 数 越 多 ， 训 练 模型 就 越 复 杂 ， 越 容易 引起 过 拟 合 
问题 。 当 两 个 染色 体 的 适应 度 函 数值 相等 时 ， 染 色 体 中 所 含 被 选中 的 特征 
个 数 少 的 优 于 所 含 被 选中 的 特征 个 数 多 的 染色 体 。 然 后 ， 确 定 遗传 算 子 产 
生 下 一 代 染 色 体 ; 利用 轮 盘 赌 的 方式 选择 适应 度 高 的 染色 体 并 复制 到 中 间 
种 群 ， 中 间 种 群 由 概率 p, =0.1 的 “一 点 交叉 ”和 概率 p,, =0.05 变异 生 
成 下 一 代 种 群 。 最 后 ， 设 定 终止 条 件 。 当 GA 的 迭代 次 数 等 于 最 大 迭代 次 
数 100 时 ， 程 序 终止 。 为 了 得 到 更 优 的 解 ， 加 入 了 一 个 基于 深度 邻 域 的 新 
改 的 方法 。 对 于 每 一 代 的 最 好 解 C, = 1…,1,…| ， 通 过 变异 每 一 位 基因 ， 
同时 其 他 基因 保持 不 变 得 到 新 的 个 体 C, = |…,0,…| ， 如 果 f(C.) >fCC,) 
时 ， 更 新 C. 为 当前 代 的 最 好 的 解 C,。SVM 是 基于 统计 学 习 理 论 的 VC 维 
(Vapnik Cherovnenkis Dimension) 理论 和 结构 风险 最 小 原理 ( Structural 
Risk Minimization Principle, SRM) [15] ; 通过 解决 一 个 线性 二 次 规划 问题 来 
求解 决策 面 ， 目 标 是 找到 训练 错误 和 泛 化 能 力 的 平衡 点 。 这 也 使 得 SVM 
方法 比 基 于 经 验 风险 最 小 化 (Empirical Risk Minimization Principle, ERM) 
的 人 工 神 经 元 网 络 等 方法 具有 更 好 的 泛 化 能 力 。SVM 将 非 线 性 可 分 样本 经 
非 线 性 变换 函数 (AK PRL) 映射 到 一 个 高 维 空间 上 ， 并 在 高 维 空间 中 构造 
最 优 分 类 超 平面 ， 这 个 超 平 面 不 但 能 将 两 类 样本 点 正确 分 开 ， 而 且 使 分 类 
间隔 最 大 。SVM 的 基本 思想 可 以 用 图 2-2 中 的 两 维 情况 说 明 ， 即 实心 点 和 
空心 点 分 别 代 表 两 类 样本 ， 为 分 类 面 ，H, 和 到 分 别 为 过 各 类 中 离 分 类 
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线 最 近 的 样本 且 平 行 于 分 类 线 的 直线 ， 它 们 之 间 的 距离 叫 作 分 类 间隔 
(margin) 。 所 谓 最 优 分 类 面 就 是 要 求 分 类 面 不 但 能 将 两 类 正确 地 分 开 〈 训 
练 错 误 率 为 0) ， 而 且 使 分 类 间隔 最 大 。 














图 2-2 ”线性 可 分 情况 下 的 最 优 分 类 务 

















SVM 的 分 类 决策 函数 为 
f(x) = sgn( È aiyik(xi,%) + b) (2-2) 
式 中 : a, 为 文 持 向 量 ; y 取 1 或 -1; 2 为 分 类 的 国 值 ; k(xi,x) =b(x,) ， 
p(x) AK PAL, 


这 里 的 分 类 器 类 型 用 的 是 C - SVM, em SUK ew iy 
斯 径 向 核 函数 和 多 层 感 知 器 核 函 数 ， 这 里 用 的 是 高 斯 径 向 核 函 数 (RBF, 
Radial Basic Function)'" ， 即 

k(x;,%;) =exp( -7 ll xx; 17) y>0 (2-3) 
式 中 : yA CH a PS tt UIE AS BY HS RK, 

SVM 分 类 器 有 一 个 缺点 ， 就 是 当 训练 样本 集中 的 两 类 样本 的 数目 相差 
很 大 时 ,分 类 结果 总 是 偏向 样本 数目 较 大 的 类 。 为 了 避免 这 一 问题 ， 本 书 
用 相同 数目 的 样本 来 训练 分 类 器 ， 而 分 类 器 的 程序 实现 是 基于 LIBSVM""7! 
完成 的 。 





























2.4 基于 SS 和 SVM 求解 模型 


在 分 散 搜 索 (SS, Scatter Search) 算法 中 ， 解 的 编码 表示 和 上 节 中 改 
HE GA 算法 的 染色 体 相 同 。 其 具体 步 又 如 下 : 
1) 随机 产生 PSize (P Size =100) 个 解 组 成 初始 种 群 P， 并 按照 解 的 
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适应 度 从 大 到 小 的 顺序 对 P 中 的 解 排 序 。 

2) 从 种 群 P 中 取 前 5 个 解 (5, =6) 加 入 参考 集 RefSet， 并 从 PP 中 去 
除 这 些 解 。 

3) 对 PP 中 的 各 个 解 分 别 求 到 参考 集 RefSet 中 解 的 最 小 距离 ， 选 择 该 
最 小 距离 最 大 的 解 加 入 参考 集 RefSet， 并 从 P 中 去 除 这 些 解 。 重 复 这 个 过 
程 5, 次 (b, =4) ， 使 得 参考 集 RefSet 的 规模 达到 6b =), + 5b,。 

4) 对 参考 集 RefSet 中 的 解 两 两 交叉 进行 组 合 ， 产 生 90 个 新 解 。 并 对 
新 得 到 的 解 进行 改进 。 

5) 用 新 得 到 的 这 些 解 更 新 种 群 P， 更 新 时 要 保留 参考 集中 的 最 好 解 。 

6) 步骤 3) ~5) 重复 执行 Maxlter 次 。 

(1) RKAS ”上面 步 又 中 需要 用 到 解 的 组 合 过 程 ， 即 利用 输入 的 两 
个 解 ， 产 生 一 个 新 的 、 具 有 这 两 解 共 同 特 征 的 解 。 在 具体 实现 时 ， 对 于 表示 
输入 解 的 2 个 二 进 制 字符 串 的 各 个 位 ， 分 3 种 情况 处 理 : 中 如 果 两 个 输入 解 
某 位 对 应 的 数值 都 为 1 ( 即 都 包含 该 特征 ) ， 则 新 解 的 相应 位 置 取 1; @) 若 都 
为 0， 则 取 0。@ 如 果 两 个 输入 解 该 位 的 数值 不 一 致 ， 则 随机 选取 0 或 1。 

(2) 解 的 距离 “SS 算法 需要 精确 度量 解 与 解 之 间 的 差异 (距离)， 这 
里 选择 文献 [11] 方法 度量 解 之 间 的 距离 ， 即 两 解 选用 的 共同 特征 个 数 与 
两 解 选择 的 特征 总 数 之 间 的 比 。 设 2 个 解 选 择 的 特征 组 合 分 别 为 X AX, , 
则 两 解 之 间 的 距离 定义 为 D(X,X,) = |X, NX, |/|X,UX,|, HP, 
|e | 表示 组 合 中 包含 的 元 素 个 数 。 

(3) 解 的 改进 SS 算法 在 迭代 过 程 中 ,需要 对 产生 的 所 有 解 进行 改 
进 。 这 里 采用 简单 的 局 部 搜索 改进 方法 ， 即 依次 对 解 的 各 位 取 反 并 计算 改 
变 后 解 的 适应 度 函 数值 ， 如 果 提 高 了 解 的 适应 度 函 数值 就 接受 改变 ， 和 否则 
该 位 改 回 原 值 。 































































































2.5 基于 改进 GA 和 SVM 模型 及 其 应 用 


本 节 以 人 体 肺 部 图 像 中 的 肺 结 节 检 测 为 例 ， 将 基于 改进 GA 和 SVM 模 
型 应 用 在 图 像 区 域 检测 中 。 图 223 是 肺 结 节 检测 的 计算 框架 ,包括 两 个 单 
元 ， 提 取 结 节 “ 候 选 点 ”和 去 假 阳 。 提 取 结 节 “ 候 选 点 ”单元 的 功能 
肺 实质 分 制 、 基 于 Hessian 矩阵 的 圆 点 滤波 器 提取 结 节 “候选 点 ”、 自 适 
应 地 选取 感 兴趣 区 域 VOL) 、 基 于 区 域 生 长 分 割 VOI 中 待 判 断 的 物体 。 
去 假 阳 单元 的 功能 有 : 提取 待 判断 物体 的 特征 、 基 于 改进 GA 和 SVM (或 
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者 基于 SS 和 SVM) 选择 最 优 的 特征 组 合 (或 者 基于 SS 和 SVM 选择 最 优 





的 特征 组 合 ) 、 用 最 优 特征 组 合 训练 出 的 SVM 分 类 器 判断 VOL 中 的 物体 是 
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2.5.1 肺 分 割 








图 2-3 肺 CAD 系统 中 的 结 节 检测 框架 








首先 ， 进 行 肺 分 割 ， 把 肺 实质 区 域 提取 出 来 ， 肺 结 节 检 测 只 需 在 肺 实 
质 区 域 中 进行 。 肺 实质 区 自动 分 割 的 方法 主要 包括 迭代 阔 值 法 、 区 域 生长 
法 、 最 优 路 径 法 和 曲率 法 。 本 节 介 绍 的 肺 实 质 分 制 方法 可 自动 、 快 速 、 准 
确 地 分 割 出 肺 实质 区 域 ， 主 要 包括 以 下 步 又: CT EGER K 
域 生长 和 连通 性 排除 肺 外 区 域 、 剔 除 气 管 、 填 充 内 部 空洞 和 剔除 肝脏 E 
右 肺 精确 分 离 和 基于 曲率 肺 区 域 修补 等 。 图 2-4a 是 肺 CT 影像 ， 图 2-4b 









































是 肺 实质 区 域 分 割 
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结果 。 





2-4 肺 CT 影像 与 分 割 肺 实质 区 域 结 
a) 肺 CT 影像 b) 分 割 肺 实质 区 
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2.5.2 自 适 应 感 兴趣 区 域 获取 


用 迭代 国 值 法 得 到 肺 实质 区 域 中 的 高 亮 区 域 互 ， 图 2-5a 表示 肺 高 亮 区 
域 的 轮廓 。 再 在 高 亮 区 域 中 用 圆 点 滤波 器 增强 具有 圆 形 结构 的 像素 点 ， 这 
些 像素 点 作为 结 节 的 “候选 点 ”集合 SCH; 图 2-5b 为 肺 结 节 “候选 点 ”。 
这 样 做 的 好 处 是 提高 速度 ， 如 果 在 肺 实质 中 所 有 像素 都 计算 Hessian 矩阵 
的 特征 值 ， 那 么 算法 的 速度 是 医生 无 法 忍受 的 。 在 5 中 进行 26 - 领域 的 连 
通 区 域 生长 得 到 多 个 候选 区 域 的 中 心 点 集合 C。 对 于 每 个 候选 区 域 提取 正 
确 特征 的 前 提 是 自 适应 地 选择 该 中 心 点 的 感 兴趣 体 的 尺寸 ， 如 果 待 判断 物 
体 的 感 兴趣 体 的 尺寸 不 对 ， 计 算出 来 的 特征 值 就 不 能 很 好 地 反映 待 判断 物 
体 的 特性 。 对 于 集合 C 中 的 元 素 ce C 其 感 兴趣 体 是 以 e 为 中 心 ， 它 的 尺 
才 从 初始 的 3 x3 x3 (像素 个 数 ) 逐步 扩大 ， 直 到 感 兴趣 体 中 高 亮 区 域 的 
个 数 占 总 感 兴趣 体 像素 个 数 的 比率 小 于 给 定 阔 值 "。 而 阔 值 + 与 “候选 点 ” 
个 数 成 正比 , > 是 0 ~0.5 之 间 小 数 。 由 于 结 节 中 是 高 亮 区 域 的 像素 个 数 
































图 2-$ 获得 自 适 应 VOI 














a) 高 亮 区 域 b) 结 节 “ 候 选 点 ” c) 结 节 VOI 
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多 ， 即 使 闵 值 + 较 大 ， 最 后 也 能 得 到 一 个 适合 的 结 节 感 兴趣 体 ， 不 至 于 因 
感 兴趣 体 太 小 而 不 能 包含 整个 结 节 。 由 于 血管 交叉 点 和 端点 处 属于 高 亮 区 
域 像 素 个 数 少 ， 即 使 闽 值 > 较 小 ， 最 后 也 能 得 到 一 个 适合 的 血管 感 兴 
体 ， 则 不 至 于 因 感 兴趣 体 太 大 而 包含 许多 与 检测 血管 无 关 的 物体 影响 判 
Wro E 2-5: 是 一 个 粘连 血管 型 结 节 三 维 可 视 化 的 结果 ， 白 色 框 为 感 兴 
体 ， 中 间 的 球状 物体 为 结 节 ， 周 边 管状 物体 为 血管 。 设 WE 是 以 e 为 中 心 
第 次 迭代 的 感 兴 趣 体 尺 寸 ，r* 是 以 c 为 中 心 第 大 次 迭代 感 兴趣 体 中 高 亮 
区 域 的 个 数 占 总 感 兴趣 体 像 素 个 数 的 比率 ， 则 有 
, #HlxeW lx eH} 
te p 





(2-4) 


2.5.3 特征 提取 


对 于 粘连 血管 型 结 节 ， 通 常 希望 只 提取 结 节 部 分 特征 ， 这 就 需要 把 血 
管 和 结 节 大 致 分 割 开 。 但 是 由 于 血管 和 结 节 的 CT 值 很 相似 ， 精 确 地 分 害 
需要 过 于 复杂 的 方法 ， 而 带 约 束 的 区 域 生长 方法 能 够 做 到 只 生长 出 结 节 的 
区 域 而 不 进入 血管 区 域 。 首 先 ， 对 $ 中 的 每 个 26 领域 的 连通 区 域 计 算 均 
Eu 和 方差 r。 然 后 ， 在 初始 候选 区 域 的 基础 上 逐渐 向 外 扩张 ， 进 行 带 约 
束 的 区 域 生长 ， 得 到 新 区 域 S。 区 域 生长 的 条 件 是 : @D 该 像素 属于 初始 
候选 区 域 像素 的 4 领域 。@ 该 像素 的 CT 值 大 于 人 -20 而 小 于 人 +2c。 并 
计算 5,., 的 均值 和 方差 。 重复 上 述 过 程 4 次 ， 最 后 得 到 的 区 域 能 被 视 作 粘 
连 血管 型 结 节 的 结 节 部 分 ， 可 以 看 到 分 割 结果 没有 渗透 到 血管 区 域 中 ， 如 
图 2-6 所 示 。 图 2-6a 为 初始 结 节 的 候选 区 域 ， 图 2-6b 为 约束 区 域 生长 后 
的 扩展 区 域 , 图 2-6c 为 VOI 中 竺 判断 物体 的 区 域 ， 即 VOI 中 的 高 亮 区 域 。 
F l 














a) b) c) 


图 2-6 特征 提取 的 3 个 区 域 
a) 初始 结 节 的 候选 区 域 b) 约束 区 域 生 长 出 的 扩展 区 域 c) VOI 中 等 判断 物体 的 区 域 
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计算 图 2-6 中 3 个 区 域 的 均值 、 方 差 、CT 最 大 值 、CT 最 小 值 和 像素 
个 数 ， 共 15 个 特征 。VOI 中 待 判断 物体 区 域 的 紧凑 度 、 球 性 、 最 小 长 短 
轴 比 率 和 对 比 度 ， 作 为 第 16 ~ 19 个 特征 。 紧 凑 度 被 定义 为 待 判断 物体 的 
内 切 球 半 径 / 外 接 球 半径 。 球 性 表示 待 判断 物体 与 其 等 体积 球 的 重 盖 度 。 
最 小 长 短 轴 比率 是 VO 中 待 判断 物体 X、Y、2Z 方向 上 长 度 任意 两 者 间 比 
例 的 最 小 值 。 上 述 3 个 特征 ， 当 待 判断 物体 是 球 时 ， 取 值 为 1， 其 他 形状 
属性 值 为 0 ~ 1 的 小 数 。 对 比 度 定义 为 带 约 束 区 域 生长 出 的 区 域 灰 度 平均 
值 与 VOI 中 竺 判断 物体 灰 度 平均 值 的 比率 。 

通过 用 服从 高 斯 分 布 的 球 来 模拟 结 节 及 用 服从 高 斯 分 布 的 圆柱 来 模拟 
血管 ， 由 于 圆柱 的 各 个 层面 像素 个 数 变化 都 比较 小 ， 而 球 在 各 个 层面 上 的 
像素 个 数 变 化 是 先 由 小 变 大 再 变 小 的 ， 则 提出 一 个 新 的 特征 作为 第 20 种 
特征 ; 计算 VOL 中 待 判 断 物体 在 XY、YZ 和 XZ 方向 上 各 层 像素 个 数 的 方 
差 ， 取 其 三 者 的 最 小 值 ， 如 图 2-7 所 示 。 
































结 节 血管 





a) 


图 2-7 各 层 像素 个 数 变 化 特征 
a) 结 节 各 层 像素 个 数 变化 情况 b) 血管 各 层 像素 个 数 变 化 情况 


























用 服从 高 斯 分 布 的 球 和 服从 高 斯 分 布 的 圆柱 的 混合 分 布 来 模拟 粘连 血 
管 型 结 节 ; 用 2 个 服从 高 斯 分 布 的 圆柱 的 混合 分 布 来 模拟 “X” 形 血管 交 
MAb; 用 3 个 服从 高 斯 分 布 的 圆柱 的 混合 分 布 来 模拟 “Y” 形 血管 交叉 
处 ， 粘 连 型 结 节 的 内 切 球 像素 个 数 /粘连 型 结 节 总 像素 个 数 的 值 会 远 远 大 
于 交叉 血管 的 内 切 球 像 素 个 数 / 交 叉 血 管 总 像素 个 数 的 值 。 当 VOL 中 待 判 
断 物 体 是 球 时 ， 取 值 为 1， 其 他 形状 特征 值 为 0 ~ 1 的 小 数 。 粘 连 型 结 节 的 
内 切 球 像素 个 数 /VOI 总 像素 个 数 的 值 会 远 远 大 于 交叉 血管 的 内 切 球 像素 
个 数 /VOI 总 像素 个 数 的 值 。 图 2-8 解释 了 上 述 2 个 新 特征 (HP, Ry 
内 切 圆 的 半径 ) ， 这 2 个 特征 为 第 21、 第 22 种 特征 。 
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图 2-8 内 切 圆 特征 
a) 粘连 血管 型 结 节 b) 立 形 血管 交叉 


2.6 实验 结果 及 分 析 


2.6.1 实验 数据 


数据 集 来 源 于 广州 医学 院 第 一 附属 医院 。 共 有 64 E CT 病例 影像 数 
据 ， 层 厚 为 2~3mm， 内 含有 结 节 共 136 个 〈 包 含 粘连 血管 型 结 节 54 个 ， 
孤立 型 结 节 82 个 ) ， 结 节 的 直径 大 于 3mm 但 小 于 3cm。 依 据 两 位 医生 共 
同 给 出 的 “ 金 标准 ”， 基 于 圆 点 滤波 器 的 “候选 点 ”提取 算法 提取 出 6253 
个 假 阳 。 分 类 器 的 训练 集 由 21 个 结 节 和 21 个 假 阳 构成 ， 而 用 于 特征 选择 
的 测试 集 由 9 个 结 节 和 439 个 假 阳 所 构成 ; 分 类 器 的 测试 集 由 全 部 的 136 
个 结 节 和 6253 个 假 阳 组 成 。 实 验 环境 和 测试 平台 都 和 第 2 章 的 相同 。 


2.6.2 参数 的 选取 


1. GA 算法 参数 的 选取 

在 进化 过 程 中 ,改进 GA 通过 训练 数据 (21 +21) 和 用 于 特征 选择 的 
测试 数据 (9 +439) 建立 的 SVM 分 类 器 来 评估 每 个 染色 体 〈 即 一 个 备 选 
的 特征 组 合 ) 的 分 类 效果 ,改进 GA 算法 的 总 迭代 数 设 为 100。 

从 表 2-1 中 可 以 看 出 ， 当 权重 参数 [IL (2-1) 释义 ] w, <0.5 时 ， 
算法 敏感 性 差 .， 这 就 使 得 训练 后 的 分 类 器 漏 检 论 高 ， 当 权重 参数 [ 见 式 
(2-1) X] w, >0.9 时 ,算法 特异 性 不 好 ， 这 就 使 得 训练 后 的 分 类 器 假 
阳 率 高 ， 所 以 改进 GA 的 评估 函数 的 权重 参数 被 设 为 w, =0.7, w, =0.3。 
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表 2-1 权重 参数 w, 和 w, 的 不 同 取 值 对 敏感 性 和 特异 性 的 影响 



































wi 0. 00 0. 10 0. 40 0. 50 0.70 0. 90 1.00 

w 1.00 0. 90 0. 60 0.50 0. 30 0. 10 0. 00 
敏感 性 0. 14 0. 93 0. 93 1.00 1.00 1.00 1.00 
特异 和 1. 00 0. 99 0. 99 0. 96 0. 96 0. 96 0. 42 
2. SS 参数 





在 产生 参考 集 时 ，N, AN, 分 别 为 评估 函数 值 好 的 解 的 个 数 和 分 散 性 
好 的 解 的 个 数 。 从 表 2-2 可 以 看 出 ，N 和 N, 取 不 同 值 时 ， 对 于 最 后 的 最 
优 解 有 一 定 影 响 。 本 节 采 用 N =6、N =4。 

表 2-2 N, FON, 的 不 同 取 值 ， 敏 感性 和 特异 性 的 变化 趋势 





























Ni 3 4 5 6 7 
N, 7 6 5 4 3 
敏感 性 1 1 1 1 1 
特异 性 0. 986 0. 986 0. 983 0. 986 0. 979 
评价 函数 值 0. 996 0. 996 0. 995 0. 996 0. 994 








2.6.3 实验 结果 及 讨论 


1. 特征 选择 算法 结果 及 讨论 
表 2-3 HERGA, WOE GA, SS 性 能 比较 。 图 2-9 是 传统 GA 和 改进 
GA 中 ， 特 征集 个 数 与 评估 函数 对 应 关系 曲线 图 。 图 2-10 是 传统 GA 和 改 





进 GA 中 ， 和 迭代 次 数 与 评价 函数 对 应 关系 曲线 图 。 











表 2-3 ”传统 GA、 改 进 GA、SS 性 能 比较 


























项 H 传统 GA 改进 GA SS 
最 优 特征 组 合 | 10100000101001101110110 | 1000000000101000110110 | 10110101001010000000110 
特征 个 数 11 7 8 
敏感 性 100. 00% 100. 00% 100. 00% 
特异 性 95. 60% 97. 50% 98. 60% 
时 间 292s 359s 151s 














从 表 2-3 中 可 以 看 出 ， 改 进 GA 比 传统 的 GA 虽然 搜索 解 的 时 间 更 长 ， 
但 搜索 的 解 更 优 (特异 性 更 好 ) 。SS 算法 比 改进 GA 搜索 解 的 时 间 短 ， 同 
时 搜索 的 解 更 优 〈 特 异性 更 好 ) 。 说 明 改 进 GA 和 SS 算法 起 到 了 作用 ， 也 
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验证 了 当 特 行 











F 不 独立 时 ， 类 别 可 分 离 判 据 不 具有 可 加 性 。 同 时 ， 最 优 解 中 


包含 用 于 区 分 血管 端点 和 孤立 型 结 节 的 新 提出 的 计算 VOI 中 竺 判断 物体 在 


R 
XY, YZ 和 XZ TA LIE RRP 22, REZA ME PA 


结 节 


包括 用 于 区 分 血管 交叉 点 和 粘连 血管 型 
数 / 待 判断 物体 总 像素 个 数 的 特 和 








FE; 还 


结 节 的 新 提出 的 计算 内 切 球 像素 个 
E ， 表 明 这 两 种 新 提出 的 特 和 





FE 的 必要 性 和 
有 效 性 。 虽 然 训练 分 类 融 的 时 间 较 长 ， 但 是 因 是 在 离线 状态 下 进行 的 ， 故 


























asg T n A.. | 
aS 
0.98 “2-8 
\ 
097 上 | 
押 
0.96 1 
a \ 
拒 0.95} ` 
lar i 
0.94 t i 4 
\ 
1 
093 上 
\ 
ag2+ oO 
0 91 L ri L L L 
6 8 10 12 14 16 18 
特征 集 个 数 
a) 
1.02 T T 
1 3 
0.98 Ao 
Fi h 
o--o 
0.96 + „E ` 
a o-- i 
E é 
= 0.94 F i 
ta Fi 1 
是 0.92r H 4 
AS 1 
09 Mg i 
! oO 
F 
ost ; 
: 
i 
ta 
0.86 区 
Oo 
0.84 1 1 : : 
6 8 10 12 14 16 18 
特征 集 个 数 





图 2-9 传统 GA 和 改进 GA H 











Ph ， 特 征集 个 数 与 评估 函数 值 对 应 关系 
a) 改进 GA b) 传统 GA 
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测试 时 间 仍 在 医生 可 以 忍受 的 范围 内 。 

从 图 2.9 中 可 以 看 出 ， 改 进 GA 在 用 较 小 特征 个 数 就 能 得 到 比 传统 CA 
还 要 好 的 效果 。 进 一 步 验证 了 ， 如 果 特 征集 合 中 的 各 个 特征 都 相互 独立 ， 
征集 合 的 个 数 越 多 ， 识 别 效果 越 好 ， 如 果 特征 集合 中 的 各 个 特征 不 相互 
独立 ,或 者 相关 性 强 ， 特 征集 合 的 个 数 就 不 一 定 越 多 越 好 。 从 图 2-10 中 
可 以 看 出 ， 改 进 CA 经 过 更 少 的 选 代 次 数 就 能 搜索 到 最 优 解 。 
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b) 


图 2-10 传统 GA 和 改进 GA F, TETRA TE Ht PRT DR HR 
a) 改进 GA b) 传统 GA 
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2. 检测 效果 及 讨论 

对 图 2-5a 的 肺 内 高 亮 区域 做 圆 点 滤波 器 操作 后 ， 产 生 了 太 多 的 候选 
点 ， 为 下 一 步 的 去 假 阳 操作 增加 了 负担 ， 也 使 得 去 假 阳 算法 的 速度 慢 。 因 
此 ， 需 要 利用 闪 值 得 掉 部 分 “候选 点 ”。 设 肺 内 高 亮 区 域 做 圆 点 滤波 器 操 
作 后 像素 最 大 的 灰 度 值 为 maxValue, WEBEN T = maxValue-k, 其 中 ke 
[0,1], Al2-lla, b 和 < 4H 24k =0.1, k=0.15 和 =0.2 时 肺 内 高 亮 
区 域 进行 圆 点 滤波 器 操作 后 的 “候选 点 ”。 可 以 看 出 ， 对 于 同一 个 提取 
“候选 点 ”算法 , k (ABR, 候选 点 就 越 少 ， 也 就 易 造成 漏 检 ; 但 假 阳 性 
率 也 少 了 ,提取 “ 候 选 点 ”算法 的 时 间 会 变 短 。 相 反 ， 如 果 〖 值 越 小 ， 候 
选 点 就 越 多 ， ee 但 假 阳性 率 又 高 了 ， 提 取 “ 候 选 
点 ”算法 的 时 间 则 会 


O E 
C) 


图 2-11 取 不 同 值 时 获得 的 “候选 点 ” 
a) k=0.1 b) k=0.15 c) k=0.2 


因为 检测 算法 只 包括 两 部 分 : 提取 结 节 “候选 点 ”算法 和 去 假 阳 算 
法 ， 所 以 漏 检 一 定 是 由 这 两 部 分 原因 造成 的 。 而 提取 “候选 点 ”导致 漏 检 
率 和 去 假 阳 导致 漏 检 率 之 和 为 1。 

“候选 点 ”提取 导致 的 漏 检 率 = “候选 点 ”提取 导致 漏 检 结 节 个 数 / 
总 的 漏 检 结 节 个 数 。 
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去 假 阳 导 致 漏 检 率 = 去 假 阳 导致 漏 检 结 节 个 数 / 总 的 漏 检 结 节 个 数 。 

忽略 “候选 点 ”因素 敏感 性 = 〈 算 法 检 出 结 节 个 数 + 提取 “候选 点 ” 
导致 漏 检 结 节 个 数 x 去 假 阳 导 致 漏 检 率 ) /总 结 节 个 数 。 该 属性 是 为 了 衡 
量 检测 算法 去 假 阳 部 分 的 效果 是 否 和 提取 结 节 “ 候 选 点 ”算法 有 关 。 

基于 表 2-4 中 的 实验 结果 ， 可 以 进行 如 下 讨论 : 

1) 由 于 该 算法 是 用 来 辅助 医生 检测 肺 结 节 的 ， 医 生还 要 对 于 算法 检 
测 出 来 的 每 个 结 节 做 出 最 后 的 结论 。 所 以 ， 医 生 更 倾向 于 漏 检 率 低 的 算 
YE, Blk =0.15 的 效果 (HE k=0.2 的 效果 ) 更 是 医生 需要 的 。 

2) 通过 分 别 分 析 天 =0.15 Fl k =0.2 的 敏感 性 、 假 阳 个 数 (E) 和 运 
行 时 间 ， 得 知 上 大 值 越 小 , “候选 点 ”的 个 数 就 越 多 ， 敏 感性 就 越 好 ， 漏 检 
率 越 低 ， 假 阳 个 数 (E) 也 就 越 多 ,运行 时 间 也 越 长 ， 反 之 亦 然 。 所 以 ， 
检测 算法 需要 取得 敏感 性 和 假 阳 个 数 (E) 的 平衡 点 。 当 然 ， 如果 值 太 
小 ， 假 阳 过 多 ， 也 会 给 医生 带 来 负担 。 

3) Kk =0.15 和 让 =0.2 的 忽略 “候选 点 ”因素 敏感 性 这 一 属性 值 基 
本 相同 可 以 看 出 ， 该 检测 算法 的 去 假 阳 部 分 算法 的 效果 与 提取 结 节 “ 候 选 
点 ”数量 这 个 参数 无 关 。 也 就 是 说 ， 该 检测 算法 的 去 假 阳 部 分 算法 的 效果 
与 提取 结 节 “ 候 选 点 ”算法 无 关 。 从 k=0.15 和 k=0.2 的 提取 “候选 点 ” 
导致 漏 检 率 和 去 假 阳 导致 漏 检 率 的 属性 值 ， 说 明 要 想 同时 提高 敏感 性 和 特 
异性 ， 需 要 从 更 好 的 提取 结 节 “候选 点 ”算法 和 去 假 阳 算法 出 发 ， 而 不 是 
单纯 的 调节 “候选 点 ”数目 这 一 参数 。 

#2-4 取 不 同类 值 时 ， 评 价 指标 值 












































项 目 k=0.2 k=0.15 

候选 点 的 个 数 / 个 4869 6253 

检 出 结 节 的 个 数 104 110 
敏感 性 76. 5% 80. 9% 
假 阳 数 / 个 〈 套 ) 0. 97 1.50 
提取 “候选 点 ”导致 漏 检 率 66. 7% 42. 8% 
去 假 阳 导致 漏 检 率 33. 3% 57. 1% 
忽略 “候选 点 ”因素 敏感 性 86. 9% 85. 2% 
平均 耗 时 /s 84. 87 121. 05 


对 于 每 个 漏 检 的 结 节 ， 分 析 其 中 的 原因 如 下 : 名 结 节 体 积 小 并 且 结 节 
上 的 像素 CT 值 暗 (在 -500Hu LAP) ， 造 成 圆 点 滤波 器 没有 滤 出 该 结 节 ， 
8 个 漏 检 结 节 属 于 这 一 种 原因 ; QP ARE, PEMA ABE 
像素 个 数 的 方差 小 ， 不 满足 第 20 种 特征 ， 这 种 原因 造成 8 个 结 节 漏 检 
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全 如 果 粘 连 血管 型 结 节 上 的 像素 CT 值 暗 (在 -500Hu 以 下 ) ， 造 成 粘连 
血管 型 结 节 初始 候选 区 域 小 ， 就 会 使 得 VOI 过 大 ， 第 21 种 特征 取 值 过 小 ， 
从 而 和 血管 交叉 点 相 混淆 ，10 个 漏 检 结 节 属 于 这 种 原因 。 

从 表 2-5 可 以 看 出 ， 飞 利 浦 产 品 的 敏感 性 虽然 高 达 100% ， 但 是 特异 
性 过 低 ， 造 成 假 阳 个 数 太 多 。GE 产品 的 优点 是 不 但 能 检测 出 实心 结 节 而 
且 还 能 检测 出 CGO 结 节 。 通 过 比较 ， 综 合 考虑 测试 数据 量 、 敏 感性 和 假 
阳 个 数 ( 套 ) ， 本 章 算法 都 达到 了 最 好 的 状态 。 

R25 ”本 章 检测 算法 与 各 产品 比较 
















































































产品 算法 数据 库 描述 fi 感 性 平均 特异 性 
waa! 25 位 病人 ，38 套 影像 ，52 个 实心 结 节 100. 0% 56. 4% 
GEUS] 50 位 病人 ，35 个 实心 结 节 和 8 个 GGO 69. 8% 每 套 影像 8. 3 个 假 阳 
GP Il) 20 位 病人 ，135 个 实心 结 节 16. 3% 每 套 影像 0. 6 个 假 阳 
本 章 算 法 64 套 病 例 ，136 个 实心 结 节 80. 9% 每 套 影像 1.5 个 假 阳 




















该 算法 不 但 可 以 检测 出 体积 较 大 的 结 节 (如 直径 在 1 ~3cm 之 间 )， 
对 于 体积 较 小 的 结 节 (如 直径 在 lem 之 下 ) 也 是 敏感 的 。 该 算法 不 但 能 
检测 孤立 型 结 节 还 能 检测 粘连 血管 型 结 证 (图 2-12 ~ 图 2-15)。 图 2-12 中 
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V 
结 节 
图 2-12 体积 较 大 结 节 的 检测 结果 
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图 2-13 ”粘连 血管 型 结 节 检 测 结果 
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图 2-14 ”孤立 型 结 节 检 测 结 
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的 结 节 为 体积 较 大 的 粘连 血管 型 结 节 ， 图 2-13 中 的 是 粘连 血管 型 结 节 ， 
图 2-14 中 为 孤立 型 结 节 ， 图 2-15 中 为 粘连 血管 型 小 结 节 。 图 2-12 ~ 
图 2-15 中 各 个 图 的 左上 、 左 下 和 右上 分 别 为 冠状 面 、 横 断面 和 矢 状 面 视 
图 ， 圈 标示 为 结 节 。 右 下 为 三 维 可 视 化 视图 结果 ， 梳 状 代表 血管 ， 团 状 代 
Ra 


























图 2-15 ”小 体积 结 节 检测 结 
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第 3 章 基于 多 尺度 和 Mean-Shift 的 
分 割 方法 研究 


3.1 多 斥 度 层次 聚 类 方法 





多 尺度 方法 是 层次 聚 类 的 代表 方法 。 多 尺度 聚 类 方法 是 一 种 模拟 
人 类 视觉 系统 的 聚 类 方法 ， 如 当 人 类 看 远 处 的 景物 时 ， 只 能 看 清楚 大 概 的 
轮廓 ， 也 就 相当 于 用 一 个 大 尺度 高 斯 模板 平滑 远 处 的 景 像 ， 把 远 处 的 影 
都 认为 属于 同一 个 聚 类 中 心 ; 当 人 类 看 近 处 的 景物 时 ， 只 能 注意 细节 ， 也 
就 相当 于 用 一 个 小 太 度 高 斯 模板 平滑 近 处 的 景 像 ， 把 近 处 的 不 同 景物 按照 
不 同 特征 归宿 到 不 同 的 聚 类 中 心 。 每 个 像素 开始 时 都 为 一 个 复 ， 当 用 不 同 
的 斥 度 平滑 图 像 时 ， 随 着 太 度 的 变 大 ， 小 的 篮 逐 渐 合 并 为 大 的 ， 直 到 整 幅 
图 像 。 多 尺度 层次 聚 类 方法 的 聚 类 过 程 如 图 3-1 所 示 。 
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图 3-1 多 尺度 层次 肾 类 树 


1. 多 尺度 理论 
多 尺度 层次 聚 类 方法 是 基于 多 尺度 理论 实现 的 。 其 图 像 和 高 斯 核 函 数 
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卷 积 的 式 为 


P(x,0) = p(x)g(x,0) = [p(x -») 





sth, p(x) 为 图 像 信号 ; g(x,0)= 1 jexp( - LS), 
p(x) F; g(x,o) (o Jim) p( F ) 


对 于 p(x) 的 极 大 值 ye R?, WKE B, = [xo E R° :limx(t,xo) =y}, # 
P, (t,x) 为 梯度 动态 系统 ，: 为 时 间 状 态 : 
dx/dt = V,P(«,o) 
Bo =K 
AP: yN B, 的 聚 类 中 心 ， 图 像 P(x) 中 的 每 个 像素 点 对 应 一 个 聚 类 中 心 。 
如 果 由 概率 密度 函数 p(x) 产 生 的 样品 集 ， 每 个 连通 的 簇 对 应 相对 高 


HE, LA PL AT AE ER, Pao) = LY 


(3-2) 

















— e|- deak TT AREER N B OLIN Parzen 估计 。 
(o /27) 20° 


根据 式 (3-2) 可 得 到 差分 方程 为 
DNOC ORTAK 
| ay 


i=l 
F exp[ = lets) -< 
fel 20° 








x(n+1) =x(n) +5; 
Cr 





(3-3) 

“eK h=o 时 ， 得 到 式 (3-4), 

2. 步骤 

1) 先 给 定 一 系列 尺度 wu，ol ，…， 且 设 oo =0 时 每 个 像素 为 一 个 
fie, BUSES UD A RRA, 

2) 由 于 尺度 o, 平 请 后 的 图 像 的 聚 类 中 心 是 由 尺度 o, ,平滑 图 像 得 到 
NM, Wo, -co = koi 假设 x(n) ARE oj 平滑 岁 像 的 聚 类 中 
心 ， 有 


a(n+1) = > o|- ee S, I? | 一 


(3-4) 
AP: x(n+1) 为 尺度 o, 平滑 岁 像 的 聚 类 中 心 ， 并 合并 具有 相同 聚 类 中 心 
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的 聚 类 子 集 为 新 的 聚 类 子 集 。 
3) 如 果 当 前 聚 类 子 集 的 个 数 等 于 或 多 于 2 个 ， 则 i=i+1， 应 重新 执 
行 第 2) 步 ; 否则 ,停止 。 





3. 定义 

聚 类 复 的 生命 期 是 指 该 徐 从 形成 到 消失 (被 合并 到 其 他 簇 中 ) 所 经 历 
的 尺度 范围 。 

4. 定理 





生命 周期 为 宽 尺 度 范围 的 聚 类 簇 比 窜 尺 度 范 围 的 聚 类 簇 更 可 能 为 真实 
物体 。 

由 此 可 知 ， 聚 类 有 效 性 的 依据 是 生命 周期 长 的 聚 类 簇 通常 对 应 于 真实 
的 物体 。 








3.2 均值 漂移 算法 


均值 漂移 算法 来 源 于 核 密度 非 参数 估计 理论 ， 它 是 基于 模式 识别 
中 的 Patten 窗 技术 。 非 参数 估计 则 是 从 样品 集 出 发 对 密度 函数 进行 估计 ， 
它 不 需要 任何 先 验 知识 ， 且 对 任何 形状 的 分 布 都 有 效 。 其 基本 思想 是 : 通 
过 反复 迭代 搜索 特征 空间 中 样品 点 最 密集 的 区 域 ， 搜 索 点 沿 着 样品 点 密度 
增加 的 方向 “漂移 ”到 局 部 密度 最 大 点 ， 即 密度 梯度 为 0 的 点 。 如 果 2 个 
羊 品 点 搜索 到 同一 个 局 部 密度 最 大 点 ， 则 说 明 该 2 个 样品 点 属于 同一 个 聚 
类 复 。 均 值 漂移 算法 中 任 一 样品 点 寻找 聚 类 中 心 的 过 程 如 图 3-2 所 示 。 
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图 3-2 均值 漂移 搜索 局 部 密度 最 大 点 的 过 程 


-32- 图 像 检测 与 分 割 方法 及 其 应 用 





1. 均值 漂移 理论 
假设 n 个 样品 点 x, i=1, =, n, x eR’, 依据 Patten 窗 理 论 ， 密 度 
f(x) = n È Rata - x) (3-5) 
式 中 : K, (x) = | 百 1 -2K( 百 -2x) 为 多 变量 核 函 数 ; HA d- d 维 的 对 称 正 
EIERE, 
H K(x) =c, k( xl W 密度 函数 式 变 为 
Fala) = > yoi ] (3-6) 


由 此 可 知 ， Tb LI a. 此 时 五 称 为 带 
WARE, WP: 61 ,为 常数 。 



















































































Saef sl 
ao Rell se] AED. 
_ ' (3-7) 
均值 漂移 向 量 为 
m, (x) = a + D, (3-8) 
Zel a 














在 上 述 式 中 ， 多 变量 核 隐 数 K(x) 必 须 满足 非 负 ， 且 积 分 为 1; 它 反 映 

的 是 样品 点 x; 与 核 中心 之 间距 离 的 度量 。 核 函数 有 多 种 ， 如 均匀 核 函数 、 

iy OT RK eK AL, ~Epanechnikov 7% pK, Biweight 核 函 数 ， 其 中 高 斯 核 隐 数 是 
最 常用 的 核 函 数 ， 其 形式 为 


















































K(x) = (2m) el -3 -Ligj i (3-9) 
当 核 函数 是 高 斯 核 函 数 时 ， 均 值 漂移 向 量 为 
Emil 
m, c(XY) = -x (3-10) 
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` x 一 Xi 
È xex | i 


a(n+1) = 一 


È ex( |: 一 
i=l 


h 











(3-11) 











2. DR 

对 于 每 个 样品 都 执行 以 下 操作 : 

1) 在 特征 空间 中 ， 以 该 样品 为 圆心 0， 半 径 为 带宽 参数 h， 得 到 初始 
搜索 区 域 圆 C。 

2) ÆR CHR (3-10) 计算 均值 漂移 向 量 m, ,(%x)， 该 向 量 总 是 指 
向 密度 增加 的 方向 。 

3) 通过 计算 圆心 0 特征 向 量 和 m, eC) 均值 漂移 向 量 之 和 得 到 新 的 
圆心 0'。 

4) 如 果 均 值 漂移 向 量 | m, ex) || <e, ERER, MAREO NZ 
样品 的 聚 类 中 心 。 

5) 与 步骤 4) 相反 ， 则 把 圆心 0' 赋 给 圆心 0， 回 到 步骤 1) 。 当 所 有 
的 样品 都 找到 聚 类 中 心 后 ， 就 会 拥有 同一 个 聚 类 中 心 的 样品 属于 一 个 聚 类 
FET 





3.3 均值 漂移 算法 带宽 参数 选取 现状 


在 Mean-Shift 算法 中 ,带宽 hh 是 一 个 重要 的 参数 ， 因 为 它 决 定 了 参与 
Mean-Shift 闪 代 的 样本 数量 ， 并 直接 影响 迷 代 速度 ， 所 以 它 对 分 割 结果 
起 着 重要 的 作用 。 如 果 和 带宽 参 数 h 选取 太 小 ， 则 偏 移 量 小 ， 算 法 收敛 速度 
很 快 ， 从 而 导致 过 分 割 现象 ， 如 果 带 宽 参 数 h 选取 大 大 ， 则 偏 移 量 大 ， 算 
法 收敛 速度 很 慢 ， 从 而 导致 欠 分 割 现象 。Mean-Shift 算法 是 基于 Patten 窗 
的 核 估计 方法 ， 核 估计 方法 既 与 同样 本 有 关 ， 叉 同 核 隐 数 K 及 徐 宽 (带宽 
BRO h 的 选取 有 关 。 在 给 定 样品 之 后 ， 一 个 核 估计 的 好 坏 ， 取 决 于 核 函 
数 及 窗 宽 的 选取 是 否 适当 。 佑 计 隐 数 在 每 个 样品 点 x, 有 一 个 “碰撞 ”， 售 
计量 是 这 些 “ 碰 撞 ” 之 和 ; 核 函 数 K 确定 了 每 个 “碰撞 ”形状 ,而 则 
决定 了 “碰撞 ”的 宽度 。 当 hh 选 得 过 大 ， 由 于 估计 函数 经 过 压缩 变换 hh 之 
后 ， 使 分 布 的 主要 部 分 的 某 些 特征 ( 如 多 峰 性 ) 被 掩盖 起 来 了 ， 估 计量 则 
有 和 较 大 偏差 ， 如 果 h 选 得 过 小 ， 这 个 估计 函数 特别 是 尾部 出 现 较 大 的 干 
扰 ， 从 而 有 增 大 方差 的 趋势 。 因 而 在 实际 使 用 核 函 数 佑 计时 ， 如 何 选 取 适 
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妆 的 宽度 是 一 项 很 细致 的 工作 。 

目前 ， 有 两 种 带宽 参数 : 固定 带宽 参数 和 可 变 带 宽 参 数 。 固 定 带宽 参 
数 是 指 所 有 样品 点 的 概率 密度 估计 都 设置 为 同一 个 带宽 参数 ， 这 种 带宽 参 
数 选 择 方法 也 叫 作 球 密度 估计 。 可 变 带 宽 参数 是 依据 数据 区 域 的 局 部 特征 
不 同 ， 为 不 同 的 点 设置 不 同 的 带宽 参数 ， 这 种 带宽 参数 选择 方法 为 可 变 带 
宽 佑 计 也 叫 作 采样 点 佑 计 ， 如 式 〈3-12) 。 下 面 的 带宽 参数 选择 定理 对 于 
这 两 种 带宽 参数 都 适用 。 基 于 统计 分 析 方 法 的 固定 最 优 带宽 参数 ， 其 评判 
规则 是 佑 计 密 度 与 真实 密度 之 间 的 误差 最 小 。 佑 计 误 差 的 度量 式 有 : 均 方 
2 [MSE， 如 式 (3-13 ) ] 、 积 分 平方 差 〈ISE) 、 积 分 均 方 差 (MISE) 、 渐 
进 积分 均 方 差 (AMISE) 等 ， 它 们 均 可 表示 为 密度 估计 函数 的 方差 与 其 偏 
差 平方 之 和 的 形式 ”i 。 通 过 分 析 发 现 ， 带 宽 参数 越 小 ， 估 计 的 偏差 越 小 ， 
而 佑 计 的 方差 则 越 大 。 因 此 应 在 两 者 之 间 进 行 折 中 ， 以 保证 估计 误差 最 
小 。 常 用 的 带宽 参数 计算 方法 有 : 快速 计算 法 、 交 又 验证 法 和 插入 规 
则 法 。 

















7 _ly_ i X-X, f 
Jamie =D AE) (3-12) 
MSE(x) =E [f(x) —f(x) ]? =Var[f(«) ] + | Bias[}(x)]} (3-13) 
Bias(x) = THp (x) fuk(u) du (3-14) 
Var(x) = ugha) [Edu (3-15) 


偏差 和 的 平方 成 正比 这 就 意味 着 小 的 带宽 参数 可 以 获得 小 的 偏差 ， 
然而 , 方差 和 成 反比 这 就 意味 着 减 小 h 将 增 大 方差 。 偏 差 表 明 佑 计量 对 
fx) 光滑 修 正 的 程度 ， 即 一 个 佑 计量 的 光滑 程度 越 高 ， 可 能 更 多 地 忽略 
f(x) 的 某 些 细节 ， 从 而 增加 随机 误差 。 一 个 核 函数 估计 的 光滑 程度 只 与 
有 关 ， 而 与 样本 集合 尺寸 无 直接 关系 。 则 使 MSE 的 值 最 小 的 有 为 最 好 
的 国定 带宽 参数 。 交 叉 验 证 、 惩 罚 函 数 、plug-in 规则 和 交叉 置信 区 间 规 则 
都 属于 统计 分 析 方 法 ， 但 它们 的 主要 缺点 是 计算 复杂 ， 因 而 不 适用 于 多 变 
量 多 模式 的 数据 处 理 。 如 plug-in 规则 中 渐进 式 给 出 了 作为 方差 和 偏差 估计 
函数 的 最 优 带宽 参数 ; 然而 ， 该 方法 依赖 于 未 知 密度 的 拉 普 拉 斯 变换 或 在 计 
算 偏差 时 引入 高 阶 导数 ， 这 都 不 适合 解决 实际 问题 。 还 有 一 种 选择 最 优 带 宽 
参数 的 方法 叫 作 交叉 置信 区 规则 后 ] ， 即 Vladimir Katkovnik 和 Iya Shmulev- 
ich 利用 假设 检验 的 自 适应 带宽 参数 选择 的 算法 。 设 H= [hy <h, < <h} 
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为 一 有 限 带宽 参数 集合 ， FEA — A AE FE BRK h MR, 确定 A(X) 的 置 
信 区 间 D(1); 依次 增 大 带宽 参数 有， 得 到 概率 密度 置信 区 间 序 列 
DG), Æ 
DQ) ND(2) NDG) Ce, H 
D) ND NDG) ND(j +1) CG 

则 由 j 值 确定 自 适应 带宽 参数 h， 并 将 该 值 作为 采样 点 概率 密度 估计 的 核 
函数 带宽 参数 。 

虽然 统计 分 析 方 法 的 可 变 带宽 可 减 小 估计 密度 与 真实 密度 函数 之 间 的 
偏差 ， 但 仍 面临 比例 参数 难以 确定 的 问题 。 选 择 使 数据 聚 类 结果 最 稳定 的 
带宽 参数 ， 是 面向 任务 的 带宽 参数 选择 方法 的 主要 思想 。 基 于 统计 分 析 方 
法 估计 带宽 参数 hh 和 实际 图 像 处 理 的 以 分 割 为 目的 的 带宽 参数 h 的 意义 是 
不 同 的 。 基 于 统计 分 析 方 法 估计 的 带宽 参数 h 是 为 了 使 佑 计量 和 实际 量 间 
的 误差 最 小 ， 以 分 割 为 目的 的 带宽 参数 (AE) h 所 反映 的 误差 程度 是 由 
实际 问题 决定 的 ， 不 是 估计 得 越 准 越 好 。 即 估计 得 准 会 使 分 割 造 成 过 分 制 
现象 ， 佑 计 得 太 不 准 了 又 会 产生 欠 分 割 现 象 。 图 像 上 以 分 割 为 目的 的 带宽 
参数 hh 是 直接 根据 样本 自 适 应 地 确定 ， 它 已 失去 了 数学 上 的 “最 佳 窗 宽 ” 
的 原意 了 。 

带宽 参数 选择 定理 天] : 假设 样本 服从 正 态 分 布 N(n,3)， 采 用 正 态 核 
函数 ky 进行 均值 漂移 ， 当 带宽 参数 h = 了 时， 均值 漂移 向 量 m(x) 的 模 取 
最 大 值 。 
































3.4 基于 多 尺度 的 带宽 参数 选取 方法 





本 章 提 出 的 结 节 精 确 分 割 算法 分 为 三 步 : 前 两 步 是 选择 初始 带宽 参数 
hoz Bes PRE RE h FH h, 代入 到 Mean-Shift 聚 类 算法 
分 割 肺 结 节 ， 得 到 精确 的 结 节 轮廓 。 首 先 利 用 2.5. 2 自 适 应 计算 感 兴趣 体 
ari eh AR 然后 用 Mean-Shift 方法 逐 层 分 割 结 节 的 区 
域 , 得 到 结 节 的 精确 边缘 。 结 合 面向 任务 的 带宽 参数 选择 方法 和 基于 多 尺 
ee ee 依据 从 粗 到 细 ( coarse-to-fine) 
准则 ， 即 从 初始 的 粗 尺度 演变 到 能 使 实际 问题 分 割 结果 最 好 的 细 尺 度 ， 提 
出 了 一 种 新 的 自 适 应 带宽 参数 选择 方法 。 
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3.5 基于 多 尺度 和 Mean-Shift 分 割 在 图 像 处 理 中 的 应 用 





本 节 用 多 尺度 和 Mean-Shift 分 割 CT 影像 上 的 肺 结 节 。 首 先 ， 利 用 区 
域 生长 方法 和 带宽 选择 定理 5 得 到 初始 的 平滑 带宽 参数 ; 然后 ， 根 据 实 
际 问题 的 特点 得 到 带宽 参数 的 最 大 值 ， 这 样 就 确定 了 带宽 参数 的 选择 范 
围 ， 再 在 这 个 范围 内 ， 依 据 多 尺度 空间 滤波 聚 类 算法 中 的 最 稳定 尺度 准 
则 ， 选 定 核 函 数 的 最 佳 带宽 参数 ; 最 后 ， 把 最 佳 带宽 参数 代入 到 二 维 
Mean-Shift 聚 类 算法 中 分 割 肺 结 节 。 基 于 频 度 和 Mean-Shift 肺 结 节 分 制程 
序 框架 如 图 3-3 所 示 。 


























结 节 检 测 




















估计 初始 带宽 参数 
优化 最 佳 带宽 参数 
Mean-Shift 


分 割 方法 


图 3-3 ”基于 多 尺度 和 Mean-Shift 肺 结 节 分 割 程 序 框架 





3.5.1 初始 区 域 的 获取 


按照 2.5.2 中 所 述 ， 得 到 肺 结 节 的 感 兴 趣 体 。 圆 点 滤波 器 提取 结 节 
REA”, FARE r 两 值 化 圆 点 滤波 器 滤波 后 的 图 像 。 图 3-4a 中 黑色 的 
点 为 结 节 “ 候 选 点 ”。 圆 点 滤波 器 滤波 后 图 像 的 直方 图 如 图 3-4b 所 示 ， 假 
设 圆 点 滤波 器 滤波 后 图 像 的 直方 图 服从 半 个 正 态 分 布 ， 自 适应 的 阔 值 r 取 
圆 点 滤波 器 滤波 后 图 像 的 方差 ， 这 种 阔 值 选取 方法 优 于 第 4 章 中 的 阔 值 
的 选取 方法 。 对 图 3-4a 中 黑色 的 “候选 点 ”做 区 域 生长 算法 ， 生 长 出 初 
始 的 区 域 ， 如 图 3-4c 所 示 。 用 本 章 的 算法 分 割 肺 结 节 的 结果 如 图 3-4d 
所 示 。 
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图 3-4 结 节 的 初始 区 域 和 本 章 的 分 割 结 果 








a) 圆 点 滤波 器 提取 “候选 点 ” b) 滤波 后 图 像 直 方 图 
c) 区 域 生长 出 初始 区 域 d) 本 章 算法 的 分 割 结果 























3.5.2 初始 带宽 参数 的 选取 


上 方 描述 的 区 域 生长 方法 还 不 能 得 到 满意 的 分 割 结果 ， 因 它 只 把 结 节 
中 较 亮 的 部 分 分 割 出 来 ， 并 没有 把 结 节 的 模糊 影 部 分 包括 进去 。 为 此 ， 需 
要 用 Mean-Shift 聚 类 分 割 方法 ， 其 中 核 函 数 的 带宽 参数 h 是 唯一 需要 确定 
的 参数 。 

根据 上 节 的 区 域 生长 的 初始 分 割 结 果 ， 确定 x 方向 的 初始 带宽 参数 
ho, y 方向 的 初始 带宽 参数 hh。， 灰 度 的 初始 带宽 参数 hoo HK RE 
长 法 分 割 出 的 结 节 中 心 高 亮度 区 域 是 服从 均值 为 ,方差 为 3 的 多 变量 
向 异性 正 态 分 布 。 根 据 带 宽 参数 选择 定理 ， 当 结 节 区 域 在 «方向 (或 y 方 
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向 或 灰 度 ) WERS dif, x 方向 (或 y 方 向 或 灰 度 ) 的 均值 漂移 带宽 参 
数 应 为 4/4， 这 是 因为 服从 正 态 分 布 WN (u, o) 物体 的 直径 等 于 4o， 因 
此 有 如 下 式 : 








X max 一 %nin = 40, = Aho 


Ymar ~ J rita =4o, =4h o (3-16) 
CT max T CT iin =4o, =4h v 
H, = diag[ hw ,hyo ,ho | (3-17) 


式 中 : xs 和 xm 为 上 节 分 割 出 的 结 节 初始 区 域 中 像素 横 坐 标的 最 大 值 和 最 
小 值 ; ys 和 yam 为 上 节 分 割 出 的 结 节 初始 区 域 中 像素 纵 坐 标的 最 大 值 和 最 
小 值 ; CT 和 CT 分 别 为 上 节 分 割 出 的 结 节 初 始 区 域 中 像素 灰 度 的 最 大 
和 最 小 CT 值 。 

当 样 本 服从 正 态 分 布 时 ， 利 用 带宽 参数 选择 定理 可 求 得 带宽 参数 h。 
当 样 本 不 服从 正 态 分 布 时 ， 利 用 基于 尺度 空间 滤波 的 多 尺度 容 类 算法 选择 
尺度 o 的 方法 来 解决 带宽 参数 选择 的 问题 。 这 是 因为 多 尺度 聚 类 算法 原理 
及 式 都 与 均值 漂移 算法 相同 。 多 扩 度 聚 类 算法 依据 聚 类 艇 的 生存 期 来 确定 
适合 实际 问题 的 尺度 参数 。 所 请 簇 的 生存 期 ， 是 指 簇 存 活 的 尺度 范围 ， 即 
从 开始 形成 该 簇 的 尺度 到 它 被 其 他 艇 合并 时 的 尺度 之 间 的 范围 。 

尺度 参数 选择 定理 : 图像 中 对 应 实际 物体 的 徐 ， 应 能 被 较 宽 的 尺 
度 范围 所 感知 。 换 言 之 ,拥有 较 长 生存 期 的 簇 比 拥 有 较 短 生存 期 的 簇 更 容 
易 被 视觉 所 感知 。 

尺度 参数 选择 定理 说 明 ， 在 基于 多 尺度 空间 滤波 的 聚 类 算法 中 ， 如 果 
在 Onin SOSO WEA, RARE (相同 或 相近 ) 的 像素 数 ， 或 者 聚 类 
簇 的 个 数 保持 不 变 ， 则 该 簇 在 尺度 ww 下 对 应 于 实际 物体 。 


3.5.3 ”最 佳 带宽 参数 的 选取 


手动 调节 参数 的 实验 表明 ， 空 间 最 佳 带宽 参数 hh, Al h, IE LET 28 
径 ， 灰 度 最 佳 带宽 hh, 正比 于 结 节 灰 度 的 变化 ,但 不 可 能 知道 这 些 先 验 知 
识 。 用 阅 值 区 域 生长 方法 分 割 结 节 时 ， 得 到 的 分 割 结果 往往 比 实际 结 节 的 
尺寸 小 ， 所 以 最 佳 带宽 参数 hh, 二 h,。 处 理 的 结 节 在 3mm ~3cm 范围 内 ， 所 
以 有 三 hh 三 h,,。 依 据 尺度 参数 选择 定理 ， 给 定 有 序 离散 带宽 尺度 集 
{hili=1,2, 1}, EP hy hh ， 用 如 下 规则 选择 一 个 最 佳 带宽 参数 
h, (AEH, hi-hi = 大 ,大 是 正 的 常数 ， 一 般 取 大 =0.1)。 

规则 1: 设 在 尺度 hh, 下， 结 节 簇 的 像素 个 数 为 C,; HRE Rh RR, 
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结 节 簇 的 像素 个 数 为 C,.,。 则 有 
minF(h,,h,,1) =C, 7C 
st hy Sh, Shyn 
规则 2: 如 果 h, 和 六 ， 都 能 使 规则 1 ABM, BFC ha,ha) = 
Fi ) ， 则 选 使 结 节 簇 像素 个 数 最 大 的 尺度 h, =max(h,,,h,) o 
规则 3: KAREN h, 时 存在 结 节 簇 中 的 像素 x*， 如 果 x% 像素 的 CT 灰 
BEAR I(x) < -7501hu。 则 最 佳 带 完 参数 < 。 
规则 1 保证 选中 的 最 佳 带 宽 参 数 h, FERRIS ASA HE; 规则 2 可 
避免 产生 局 部 极 小 值 ; 规则 3 可 避免 选 i 党 肺 组 织 聚 到 一 起 的 
尺度 ， 避 免 了 欠 分 割 现象 。 将 选中 的 最 佳 带 宽 参 数 h, 代入 Mean-Shift 聚 类 
算法 , 便 可 得 到 结 节 簇 。 
3.5.4 Mean-Shift 二 维 分 割 结 节 过 程 
节 感 兴趣 区 域 中 的 像素 映射 到 空间 位 置 特 征 和 灰 度 特征 的 联合 特征 
空间 中 , X, = (x,y, eR 并且" =1(x,,y7,), 了 表示 灰 度 ; 然后 用 


Mean-Shift 对 感 兴趣 区 域 逐 进行 从 类 分 割 。 基于 高 斯 核 的 密度 函数 估计 为 
1 


P 


(3-18) 



































RD =e eee > epf- FP (X,X,,H)| (3-19) 
H = diag{ h* ,h’ ,h’] (3-20) 
D(X,X,,H) =(X -X,)"H"'"(X -X,) (3-21) 


式 中 ; HA3 x3 的 对 角 带 宽 参 数 矩 阵 ，D 为 两 个 样品 在 特征 空间 上 的 马 
氏 距 离 。 
对 7(X) 求 密度 梯度 : 


a H” (X, - X) 
SO rank HI” exp - FD°(X,X,H) | (3-22) 





È xep- FPL, H) | 
m(X) = =X (3-23) 
È epf- 了 Oo 
Mean-Shift 向 量 总 是 指向 密度 增长 最 快 的 方向 。Mean-Shift EM Wik 


代 下 面 中 、@ 两 步 ， 保 证 收敛 到 概率 密度 估计 函数 梯度 为 0 的 点 。 
D 计算 Mean-Shift 向 量 m(X') ; 
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Q PRA, BAP OWE X, Xt =X 4+m(X'), 

该 算法 迭代 地 分 割 肺 结 节 边 缘 的 步骤 如 下 : 

1) 初始 化 (t=0): 

D 求 结 节 感 兴趣 区 域 ， 并 区 域 生长 出 初始 区 域 R ; 

© 用 式 (3-16) 计算 初始 带宽 参数 hy 

2) 迭代 (t=t+1): 

D Fh, -h,_., =kh, ERTES h; 

@) Ki h, 代入 到 均值 漂移 算法 中 分 市 感 兴趣 区 域 ， 得 到 区 域 R,; 

© 如 果 区 域 R 满足 3.5.3 所 有 规则 ， 则 最 优 带宽 参数 h, =h, A 
步骤 3); 否则 ， 再 回 到 步骤 2)。 

3) 最 后 ， 将 最 优 带 宽 参 数 h, 代入 到 基于 均值 漂移 算法 中 分 割 感 兴趣 
区 域 ,得 到 精确 的 肺 结 节 边缘 。 











3.6 实验 结果 及 分 析 


3.6.1 实验 方案 


本 章 的 实验 环境 和 实验 平台 都 与 第 3 章 相同 。 本 章 共 设计 了 3 BK 
验 ， 具 体 方案 如 下 。 

实验 一 : 设计 人 造 数据 实验 ， 目 的 是 测试 新 分 制 算 法 的 抗 噪声 的 
实验 二 : 设计 肺 CT 数据 实验 ， 目 的 是 测试 不 同 的 带宽 参数 ， 聚 类 分 
割 结果 是 不 同 的 。 用 不 同 的 带宽 参数 的 Mean-Shift 算法 分 割 肺 结 节 ROT, 
在 空间 位 置 关系 — 灰 度 特征 直方 图 上 观察 分 割 结果 。 

实验 三 : 用 新 分 割 算 法 分 割 实验 数据 中 的 肺 结 节 ， 分 割 结果 与 “ 金 
标准 ” 相 比 较 。 带 宽 参 数 选择 范围 与 文献 [27] 的 带宽 参数 选择 范围 相 
比较 。 
3. 6.2 人 造 实验 数据 及 结果 

为 了 验证 本 章 分 割 算法 抵抗 噪声 的 能 力 ， 构 造 二 维 图 像 ， 图 像 是 灰 度 
值 服从 各 向 异性 的 高 斯 分 布 的 椭圆 。 通 过 不 断 地 给 人 造 的 图 像 加 入 噪声 ， 
并 用 本 章 算法 分 割 被 噪声 污染 后 的 图 像 ， 最 后 分 割 结果 和 “ 金 标 准 ” 做 比 
BE, 来 测试 本 章 算法 对 噪声 的 敏感 程度 。 图 3-5a 是 灰 度 服从 正 态 分 布 
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Pe =2 
N(u,2) HEN, JEP u = [0,0] HNP, Irs x =|” 5 pat 
圆 的 跨度 。 分 割 结果 的 “人 金 标准 ”是 图 3-$a 的 50% 置信 椭圆 ， 如 图 3-5b 























所 示 。 用 不 同 太 度 的 高 斯 噪声 (高 斯 模板 ) 和 图 3-5a 进行 卷 积 处 理 已 达 
到 往 图 像 上 不 断 加 入 高 斯 噪声 的 效果 。 用 本 章 算 法 分 割 受 噪 声 污染 的 图 
像 ， 可 以 计算 分 割 结果 的 错误 率 。 图 3-5c 所 示 的 横 轴 是 噪 









































声 的 级 数 ， 也 























就 是 高 斯 尺度 参数 的 取 值 ， 纵 坐标 是 分 割 的 错误 率 。 当 高 斯 尺度 参数 为 











[0,0.5] 























时 ， 分 割 加 入 噪声 后 的 图 像 错 误 率 为 0， 即 分 割 结 

















高 斯 尺度 


的 增加 ， 对 图 像 的 噪声 污染 程度 也 增加 了 ， 分 割 错 











大 。 从 图 3-5c 中 可 以 看 出 ， 当 高 斯 尺度 从 0 增加 到 0.5 AY, 





正确 ， 








说 明 该 算法 具有 一 定 的 抗 噪声 能 力 。 





a) b) 


分 割 的 错误 率 
© 
A 





0 05 1 is 2 
噪声 的 级 数 
c) 
图 35 ”人 造 数据 实验 结 


a) 2D 高 斯 人 造 椭圆 b) 分 割 结果 的 “ 金 标 准 ” 
c) 噪声 和 分 割 错 误 率 对 应 关系 

















REH, ME 
误 率 也 越 来 越 
分 割 结果 仍然 
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3.6.3 肺 CT 影像 实验 数据 


评估 实验 的 CT 数据 来 源 于 各 向 同性 的 LIDC 数据 库 和 各 向 异性 的 广 
州 医学 院 第 一 附属 医院 数据 库 。 共 18 个 病例 ，36 个 肺 结 节 ， 共 95 幅 含有 
结 节 的 CT 影像 ， 其 中 : 有 7 幅 含 有 毛 玻 璃 型 结 节 的 影像 ，21 幅 含 有 粘连 
肺 壁 型 结 节 的 影像 ，27 幅 含 有 粘连 血管 型 结 节 的 影像 ，40 幅 含 有 扳 立 型 
结 节 的 影像 。 


3.6.4 实验 结果 与 讨论 


图 3-6 是 分 析 带 宽 参 数 对 Mean-Shift 分 制 结果 的 影响 。 可 以 看 出 ， 当 
用 不 适合 的 带宽 参数 时 ,分 割 结果 会 造成 过 分 割 和 欠 分 割 现象 。 图 3-6a 
是 图 3-6c 中 的 结 节 ,图 3-6b 是 图 3-6a WARE -空间 位 置 关系 二 维 直 方 
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图 3-6 不 同 带宽 参数 聚 类 结果 二 维 直 方 图 
a) VOI 中 的 结 节 b) Ala 的 二 维 直方 图 
c) 初始 带宽 参数 分 制 结果 直方 图 d) 最 佳 带宽 参数 分 割 结果 直 方 图 
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图 ,x 和 yy 方向 表示 空间 位 置 特征 ， 而 z 方 向 表示 灰 度 特征 。 图 3-6c 是 带 
宽 参数 为 初始 带宽 参数 时 的 分 割 结果 ， 可 以 看 出 没有 把 结 节 上 的 像素 全 部 
聚 成 一 个 复 ， 从 而 产生 了 过 分 审 。 图 3-6d 是 带宽 参数 为 最 佳 带宽 参数 时 
的 分 割 结 果 ， 即 把 结 节 上 的 像素 全 部 聚 成 一 个 复 ， 显 示 分 制 结果 正确 。 

区 域 生 长 方法 的 分 割 结果 如 图 3-4c 所 示 ， 用 本 章 算法 的 分 割 结果 如 
图 3-4d 所 示 。 比 较 图 3-4c 和 d 的 结果 可 以 看 出 ， 传 统 区 域 生长 方法 只 能 
分 割 结 节 中 心 的 高 亮 区 域 ， 却 不 能 包含 结 节 周 围 的 模糊 影 ; 而 本 章 算法 则 
能 精确 分 割 结 节 边缘 ， 且 其 结果 符合 医生 手工 勾画 出 的 边缘 。 

图 3-7 是 对 含有 不 同类 型 结 节 的 分 割 结果 ， 表 明 该 算法 对 毛 玻 璃 型 结 
节 、 孤 立 型 结 节 (包括 直径 小 于 5mm 的 结 节 ) 、 粘 连 肺 壁 型 结 节 、 粘 连 血 
管 型 结 节 都 是 有 效 的 。 但 是 对 于 空洞 型 结 节 该 算法 是 无 效 的 。 
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在 相同 的 迭代 步 长 k=0.1 下 ， 本 章 算 法 分 割 如 图 3-7 中 的 4 个 结 节 。 
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参数 选择 范围 与 文献 [27] 中 算法 在 分 割 同样 的 这 4 个 结 节 时 ， 
数 选择 范围 的 对 比 见 表 3-1。 可 以 看 出 ， 本 章 算法 大 大 缩小 了 带 
选择 范围 ， 也 即 减 少 了 迭代 次 数 。 如 对 于 图 3-7a 中 的 孤立 型 小 
2a, SCHR [27] 算法 的 最 佳 带 宽 参 数 选 择 范围 为 0.5 ~ 13 ， 而 本 章 算法 
的 最 佳 带宽 参数 选择 范围 为 1. 500 ~4. 025， 缩 小 了 大 约 四 分 之 一 。 


表 3-1 带宽 选择 范围 的 对 比 


















































图 号 文献 [27] 算法 本 章 算法 
R| 3-7a 0.5~13 1. 500 ~4. 025 
R| 3-7b 0.5 ~13 1. 500 ~ 10. 500 
RI 3-7 0.5 ~13 3. 000 ~9. 750 
| 3-7d 0.5 ~13 4. 250 ~6. 825 

















第 4 章 基于 EM 和 Mean-Shift 
的 分 割 方 法 研究 


4.1 期 望 最 大 算法 


4.1.1 算法 原理 

(BSE REAR JF SRR HE CO =1,2,…,m) 的 概率 为 P(C,)， 对 应 各 个 艇 
的 条 件 概 率 密度 函数 为 p(x1C;,0,)， 其 中 b 为 未 知 参数 向 量 。 这 样 由 m 
类 样本 组 成 的 混合 密度 函数 可 定义 为 








p(x10) = È pll CG,0)P(C) (4-1) 


式 中 : > PCC) = 1; 混合 概率 密度 的 参数 0 = (6, ,0,,… ,0,)。 
7=1 


最 大 似 然 概率 (ML, Maximum Likelihood) [26] 是 参数 估计 的 一 种 重要 
的 方法 ， 即 对 n MER X = 4x, ,x,,…,x,| ， 求 出 适当 的 参数 ， 使 联合 概率 
密度 函数 式 (4-1) 取得 最 大 值 。 





pdam 1 0) = Dp! 0) (4-2) 
其 对 数 的 形式 为 


KODO = WT Tr(a! 0) = Yie( LPC rem!) 043) 





通过 (91(0) )/(00,) =0 求 得 参数 0 = arg max] [pC 0) 
遗憾 的 是 ， 最 天 似 然 估计 法 在 天 多 数 情况 下 不 能 用 解析 法 求解 ， 因 而 
需要 迭代 的 方法 寻找 优 解 。 期 望 最 大 算法 就 是 其 中 最 常用 的 一 种 方法 。 
用 EM 处 理 “ 不 完整 ”的 样本 ， 即 把 样本 点 分 为 两 部 分 ,x = aa ho 
x 表示 样本 ，x,”= (x a) 表示 丢 失 的 数据 。 岂 是 1 还 是 0， 是 根据 x 
否 属于 类 i 有关!'” 。 完 整数 据 概率 对 数 形式 为 
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(0) = X È xsl P(C) p(s | 0,)] (4-4) 


j=l i=l 


用 EM 估计 混合 概率 密度 的 参数 时 ， 需 要 首先 给 出 参数 的 初始 解 0 ; 
再 从 初始 解 8" 开始， 迭代 地 得 到 解 9” ，…，6” 。 在 每 步 迭 代 中 ， 似 然 
函数 单调 增加 。 


4.1.2 算法 步骤 


基于 以 上 算法 原理 和 公式 ， 具 体 步 又 如 下 : 

1) 给 定 初始 解 0 , + =0。 

2) MFral, 2, =, E-W, FEMUR PEAS ASAT O° AAEM ARE FP, 
计算 完全 数据 的 对 数 似 然 函 数 期 望 值 ， 根 据 式 (4-4) 有 

0(0,0° ) =E[lgp(x’,x"|10) |x ,0"] (4-5) 

3) M- 步 :选择 一 个 新 的 参数 使 得 0- 函 数 最 大 ，0 0) = arg maxsQ(9， 
0 )。 

4) t=t+1, 重复 步骤 2)、3) ， 直 到 条 件 收敛 。 








4.2 基于 EM 和 Mean-Shift 在 图 像 分 割 中 的 应 用 


4.2.1 基于 EM 的 带宽 选取 方法 流程 


本 节 将 EM 和 Mean-Shift 算法 应 用 到 3DCT 粘连 血管 型 肺 结 节 的 分 割 
中 。 根 据 文献 [28] 的 算法 的 步骤 可 以 推断 出 结 节 是 服从 高 斯 分 布 的 球 ; 
血管 是 服从 高 斯 分 布 的 圆柱 ; X 形 血 管 是 由 2 个 服从 高 斯 分 布 的 圆柱 组 
成 ; Y 形 血管 是 由 3 个 服从 高 斯 分 布 的 圆柱 组 成 。 由 此 也 可 以 推断 出 粘连 
血管 型 结 节 是 由 服从 高 斯 分 布 的 球 和 服从 高 斯 分 布 的 圆柱 组 成 。 

血管 里 充满 着 向 同一 方向 流动 的 血液 ， 血 管 呈 管 状 结构 ， 通 常 把 血管 
模拟 为 服从 高 斯 分 布 的 圆柱 ， 这 就 使 得 血管 上 像素 梯度 的 垂直 方向 基本 一 
致 ， 血 管 上 像素 梯度 的 垂直 方向 集合 服从 正 态 分 布 ， 而 结 节 里 充满 着 癌 细 
胞 ， 结 节 是 球状 结构 ， 通 常 把 结 节 模 拟 为 服从 高 斯 分 布 的 球 ， 这 就 使 得 结 
节 上 像素 梯度 的 垂直 方向 指向 四 面 八方 ， 结 节 上 像素 梯度 的 垂直 方向 集合 
服从 均匀 分 布 ， 如 图 4-1 所 示 。 

首先 利用 关系 矩阵 求 得 梯度 的 法 向 量 方向 ， 即 流向 特征 ， 然 后 ， 假 设 
粘连 血管 型 结 节 像素 梯度 的 垂直 方向 服从 正 态 分 布 和 均匀 分 布 的 混合 概率 
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a) b) 
图 4-1 结 节 和 血管 上 各 像素 的 梯度 方向 
a) 结 节 梯 度 方向 b) 血管 梯度 方向 








分 布 ， 采 用 期 望 最 大 算法 拟 合 出 参数 和 权重 。 均 匀 分 布 的 权重 可 以 被 看 作 
是 结 节 上 像素 个 数 占 粘 连 血管 型 结 节 总 像素 个 数 的 比重 。 最 后 ， 根 据 均 义 
分 布 的 权重 求 得 核 函 数 的 带宽 参数 ， 代 入 均值 漂移 ( Mean-Shift) 聚 类 分 
割 出 精确 的 肺 结 节 边缘 。 人 造 数据 和 肺 结 节 真实 数据 的 实验 结果 表明 ， 该 
方法 快速 而 有 效 地 解决 了 粘连 血管 型 结 节 的 分 割 问题 。 

肺 结 节 分 割 框架 如 图 4-2 所 示 。 首 先 确 定 粘连 血管 型 结 节 的 感 兴趣 区 
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提取 流向 特征 


建立 粘连 血管 型 结 节 模 型 


输出 带宽 参数 


Mean-Shift 三 维 聚 类 
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域 (VOI) 。 然 后 对 感 兴趣 区 域 (VOL) 进行 预 处 理 ， 得 到 粘连 血管 型 结 节 
的 前 景区 ， 然 后 提出 前 景区 的 流向 特征 ， 最 后 建立 粘连 血管 型 结 节 模 型 。 
估计 模型 参数 ， 根 据 估计 出 的 模型 参数 输出 带宽 参数 ， 带 入 到 均值 漂移 算 
法 中 进行 三 维 粘连 血管 型 结 节 分 割 。 


4.2.2 预 处 理 


对 于 图 4-3 的 粘连 血管 型 结 节 ， 做 如 下 操作 。 首 先 ， 用 2.5.2 中 的 方 
法 自 适 应 地 得 到 结 节 感 兴趣 体 ; 然后 ， 用 和 迭代 阔 值 法 得 到 的 阔 值 进行 
26- 领 域 区 域 生 长 的 方法 分 割 粘 连 血管 型 结 节 的 前 景区 〈 即 包含 结 节 和 血 
管 部 分 ) ; 最 后 ， 用 最 大 密度 投影 (MIP, Maximum Intensity Project) 向 横 
断面 方向 做 投影 ， 投 影 到 二 维 平面 上 ， 如 图 4-4 所 示 。 最 大 密度 投影 是 指 
对 于 三 维 感 兴趣 体 按 某 个 方向 CX 轴 方 向 、Y 轴 方 向 或 Z 轴 方 向 ) 投影 到 
一 个 二 维 平面 上 。 二 维 图 像 的 CT 值 为 各 个 层 对 应 位 置 像 素 CT 值 的 最 大 值 。 



































图 4-3 ”一 个 二 维 血管 粘连 型 结 节 








图 4-4 最 大 密度 网 4-3 的 结 节 
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4.2.3 基于 关系 矩阵 的 流向 特征 提取 


文献 [28] 中 给 出 了 三 维 的 关系 矩阵 ， 这 里 导出 二 维 的 关系 矩阵 。 在 
MIP 图 像 上 (图 4-4) 以 像素 j 为 中 心 的 3 x3 窗 中 的 二 维 高 斯 梯度 向 量 
6 = |gi1 1， 其 中 n=3 ”， 由 于 ROI 尺寸 小 ， 高 斯 尺度 也 应 该 小 。 垂 直 于 
6 的 向 量 v 被 定义 为 投影 到 中 所 有 向 量 的 平方 和 最 小 E(v) = 




















Z 2 1 3 
+> ((g;) y =v" GG"v, 其 中 ， G= lag GG' 是 2 x2 关系 
vt n 


和 矩阵。A, A, AGG" 的 特征 值 且 A, <A,, SA, 对 应 的 特征 向 量 是 e ， 则 
像素 i 的 梯度 的 法 向 量 为 

v,= VA, xe, (4-6) 
RP: v ALAA ta], WME 4-5 所 示 。 
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图 4-5 流向 特征 图 





图 4-5 是 图 4-4 的 流向 特征 图 。 虽 然 特征 向 量 有 唯一 的 方向 ， 但 是 特 
征 向 量 并 不 是 唯一 的 。 为 了 让 向 量 充满 [0,27] 区 间 ， 需 要 重新 定义 特 
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征 向 量 的 象限 : 





Ifg,(% 5%) E I Thenv,(x),¥%) € I 
Ifg,(% 5%) € I Thenv,(x,,¥) € M (4-7) 
Ifg,(% 5%) € MThenv;cxo, yo) € V 
Ifg,(% 5%) € WThenv,(x%),7¥)) € I 
为 了 保证 ROI 区 域内 血管 流向 一 致 ， 需 要 把 一 些 向 量 的 方向 转换 到 其 
相反 的 方向 。 在 指向 同一 方向 最 多 的 向 量 中 任 取 一 个 向 量 Arbv， 如 果 它 和 
其 他 向 量 ”之 间 的 夹 角 在 人 邢 ,] 范 围 内 ， 把 向 量 ” 转 到 相反 的 方向 -， 
如 图 4-6 所 示 。 
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图 4-6 一 致 化 流向 特征 图 


v, + Arby 
ie —y, "Ln (4-8) 
v, ”其 他 
对 向 量 的 方向 角度 做 直方 图 。 如 果 直 方 图 不 连续 ， 把 部 分 角度 加 2T 
后 使 得 直方 图 连续 并 所 有 角度 在 [0,7] 范围 内 ， 如 图 4-7 所 示 。 
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图 4-7 图 4-6 的 直方 图 





4.2.4 ”基于 EM 估计 粘连 血管 型 结 节 模 型 参数 


假设 ， 结 节 是 服从 高 斯 分 布 的 球状 物 ， 血 管 是 服从 高 斯 分 布 的 管状 物 ， 

则 结 节 梯 度 的 法 向 量 方向 指向 四 面 八方 的 ， 服 从 均匀 分 布 50, 100], m 

管 梯度 的 法 向 量 方向 是 沿 血管 流 向 指向 的 ， 服 从 正 态 分 布 ， 如 图 4-7 所 
示 。 因 此 ， 建 立 粘 连 血管 型 结 节 的 模型 如 下 : 

E 1 _ (x-u)*) os 

Ka) =o, ep -Eaa (4-9) 

RP: o, 和 w, 为 权重 并 有 o, +o, =1; y Mo 分别 为 正 态 分 布 的 均值 和 

方差 ; x 为 前 景区 (二 维 平面 ) 各 个 像素 流向 特征 角度 ; f(x) 为 角度 的 分 



































根据 第 5 章 介绍 的 带宽 参数 选择 的 方法 有 : 以 c AP, 尺寸 为 
rxy! rt rz 的 3DVOI 数 据 的 特征 空间 XX, = (x,,y,,z,,r;) ER’, 其 中 = 
Tx;,y;,2;) 5 Xs Vix z 为 空间 位 置 特 征 ,， r, 为 灰 度 特征 。 已 知 MIP 的 前 景 
区 中 的 像素 个 数 NV,, ， 根 据 均匀 分 布 的 权重 w, ， 可 知 MIP 的 前 景区 中 结 节 
TBAT HY RA A BUEN, w,。 结 节 在 x、y 方向 上 的 跨度 是 Sy = Round 
(Nu %;)， 在 z 方 向 上 的 跨度 是 S, = Round (S,, wj/zwmm)，Round 为 
BERE, £un I CT 图像 x、y 方向 的 分 辩 率 ，zw 为 z 方 向 的 分 辩 率 ,在 
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方向 上 的 跨度 是 $ = CT - CT ，，C7T，，67T 分 别 为 N,, 中 像素 的 最 大 
和 最 小 CT 值 。 假 设 结 节 服 从 各 向 异性 的 正 态 分 布 ,，x 方向 的 带宽 参数 h., 
y 方 向 的 带宽 参数 六 ，z 方向 的 带宽 参数 h, 和 灰 度 的 带宽 参数 h,， 根 据 带 
宽 参 数 选择 定理 ， 当 结 节 区 域 在 x 方向 (或 y、z 方 向 或 灰 度 ) 的 直径 为 d 
时 , x 方向 (或 y、z 方 向 或 灰 度 ) 的 均值 漂移 带宽 参数 应 为 w4， 这 是 因 





为 服从 正 态 分 布 N(p,o ) 物 体 的 直径 等 于 A40, 
h, =S,,/4 
h,=S,/4 
7 a (4-10) 
h, =S/4 
h, =S,/4 





4.3 Mean-Shift 三 维 分 割 结 节 过 程 


对 于 rxy" t rxy" : rz" 的 三 维 感 兴趣 体 数据 集 X; eR’, tly 14; N us 
MPERA X AE BK) EZE E AER E E RRO 














Nou 





7 E 1 x o1 j y 
[= ae p[- SPX] 11) 
D(X, X,,H) =(X =X)" .H '(X-X.) (4-12) 
H = diag(h, ,h, ,h,,h,) (4-13) 








寺 征 空间 中 ， 密 度 的 极 大 值 点 为 聚 类 中 心 点 。 也 就 是 ， 密 度 函 数 梯度 
等 于 零 的 点 ， 即 

1 人 再” (X, -X) 
N (27)° IHI 


wu 


Vx) = exp| -D(X,X,H) | (4-14) 





Nuu 
> Xexp[- 5-0°(4,%,,H) | 
m(X) = 5 


a T 
È exp -二 PC 
为 聚 类 中 心 过 程 中 的 偏 移 向 量 。 
Mean-Shift 向 量 总 是 指向 密度 增长 最 快 的 方向 。Mean-Shift st Fei wh ik 
代 下 面 1)、2) PAA, PRUE TUBERS AD 
1) 计算 Mean-Shift 向 量 m(X) 。 
2) 平移 窗口 ， 窗 口中 心 位 置 移 到 X, X =X +m(X')。 结 束 条 件 
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BX =X, RNBAAARAN ES, FIRER, X 是 样本 在 第 上 次 迭 
代 中 的 坐标 和 CT 灰 度 值 ，X ”是 样本 在 第 t+1 次 迭代 中 的 坐标 和 CT 灰 
度 值 。 对 于 每 个 样品 点 ， 运 行 上 述 过 程 ， 收 敛 到 相同 容 类 中 心 的 点 的 集合 
被 认为 属于 同一 个 聚 类 复 ， 从 而 得 到 精确 的 结 节 轮 廓 线 。 


4.4 结果 及 分 析 


4.4.1 实验 方案 


本 章 的 实验 环境 和 实验 平台 都 与 第 3 章 相 同 。 本 章 共 设计 了 2 BK 
验 ， 具 体 方案 如 下 。 

实验 一 : 设计 人 造 数 据 实验 ， 目 的 是 测试 流向 特征 是 否 能 区 分 孤立 型 
结 节 、 粘 连 血管 型 结 节 、 血 管 端点 、X 形 血管 交叉 点 和 工 形 血管 交叉 点 。 

实验 二 : 设计 肺 CT 数据 实验 ， 用 新 分 割 算法 分 割 ， 并 三 维 可 视 化 实 
验 数据 中 的 肺 结 节 。 分 割 结果 与 “ 金 标准 ” 相 比 较 ; 带宽 选择 范围 与 文献 
[27] 和 第 5 章 的 带宽 参数 选择 范围 相 比较 。 


4.4.2 人造 实验 数据 及 结果 


1 个 半径 为 6 (像素 ) 、 尺 度 为 1 的 三 维 服从 高 斯 分 布 的 球 与 1 个 半径 
为 2 (像素 ) 、 高 度 为 21 (像素 )、 尺 度 为 1 的 三 维 服从 高 斯 分 布 的 圆柱 相 
用 加 模拟 粘连 血管 型 结 节 。 把 三 维 物体 最 大 密度 投影 ， 如 岁 4-8a 所 示 ; 
求 出 图 4-8a 的 流向 特征 向 量 ， 如 图 4-8b 所 示 ; 流向 特征 向 量 方向 直方 图 
如 图 4-8c 所 示 ， 横 轴 为 流向 特征 向 量 方向 角度 ， 范 围 为 [0,T] ， 纵 轴 为 
像素 个 数 。 结 果 显 示 直 方 图 服从 均匀 分 布 和 高 斯 分 布 的 又 加 。 

1 个 三 维 服从 高 斯 分 布 的 球 模拟 孤立 型 肺 结 节 最 大 密度 投影 后 的 流向 
村 征 向 量 如 图 4-9al 所 示 ， 流 向 特征 向 量 方向 直方 图 如 图 4-9a2 所 示 。 可 
见 直 方 图 服从 均匀 分 布 。 

1 个 三 维 服从 高 斯 分 布 斜 圆 柱 模拟 血管 最 大 密度 投影 后 的 流向 特征 向 
量 如 图 4-9bl 所 示 ， 流 向 特征 向 量 方向 直方 图 如 图 4-9b2 所 示 。 可 见 直 方 
图 服从 高 斯 分 布 。 

2 个 三 维 服 从 高 斯 分 布 的 圆柱 相 徐 加 模拟 X 形 血 管 最 大 密度 投影 的 流 
向 特征 向 量 如 图 4-10al 所 示 ， 流 向 特征 向 量 方向 直方 图 如 图 4-10a2 所 示 。 
直方 图 服从 2 个 高 斯 分 布 相 炙 加 。 
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a) 人 造 粘 连 血管 型 结 节 MIP 结果 b) 流向 特征 向 量 c) 流向 特征 向 量 方向 角度 直方 图 


mR 


























a N 
pyra 
AILLI om, 
AAA 
MEET A EA 
SISSIES 
Po 
parr 
AIAAUAUIAAIYAT 
AAAIAAIAAIAIAAT 
YUHIYUIGGYYYT 
Po 
fons VIAAAAVAAIAAIAY 
es AAAAAAAAAUAUAA 
Fon Rt 
f AA 
a 7A7AAAAAAAAAT 
Na ed AAAAIAAIAIIAIT 
es ed SII 
MY reer 了 
Sas ceri RRA AARR ER RARE 
AVIUIIIIYIIIYZY 
了 AAAAAUIAAAAZT 
ia 
了 zz 7 
NW 
Ne er 


1200 
1000 











a2) b2) 

图 4-9 ”人造 孤立 型 结 节 、 血 管 流向 特征 向 量 及 其 方向 直方 图 
al) 人 造 抓 立 型 结 节 流 向 特征 向 量 ， bl) 人 造血 管 流向 特征 向 量 
Ja 量 


a2) 人 造 扳 立 型 结 节 流 向 特征 向 量 方向 直方 图 b2) 人 造血 管 流向 特征 向 
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3 个 三 维 服 从 高 斯 分 布 的 圆柱 模拟 Y 形 血 管 最 大 密度 投影 后 如 
图 4-10bl 所 示 ; 流向 特征 向 量 方向 直方 图 如 图 4-10b2 所 示 。 可 见 直方 图 
服从 3 个 高 斯 分 布 禾 加 。 
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图 4-10 人 造 X 形 和 YY 形 血 管 流 向 特征 向 量 及 其 方向 直方 图 
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a2) 人 造 X 形 血管 流向 特征 向 量 方向 直方 图 b2) 人 造 Y 形 血管 流向 特征 向 量 方向 直方 图 





4.4.3 肺 CT 影像 实验 数据 


评估 实验 数据 全 部 来 源 于 各 向 同性 的 LIDC 数据 库 和 各 向 异性 的 广州 
医学 院 第 一 附属 医院 数据 库 。 其 中 ，LIDC 数据 库 中 包含 6 个 病例 中 的 6 
个 粘连 血管 型 结 节 ， 数 据 库 中 包含 专家 分 割 结 节 的 “ 金 标准 ”; 广州 医学 
院 第 一 附属 医院 数据 库 中 包含 10 个 病例 中 的 10 个 粘连 血管 型 结 节 ， 数 据 
库 中 包含 医生 画 定 的 结 节 分 割 结 
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4.4.4 实验 结果 与 讨论 


采用 本 章 算法 分 割 这 16 个 粘连 血管 型 结 节 ， 把 算法 分 割 的 结果 和 医 
生 的 “ 金 标 准 ” 相 比较 ， 其 结果 全 部 符合 医生 手工 勾画 出 的 边缘 。 图 4-11 
是 对 图 42 中 的 粘连 血管 型 结 节 的 各 层 分 割 结 果 。 与 区 域 生 长 算法 相 比 
较 ， 可 以 看 出 ， 本 章 算 法 能 把 粘连 血管 型 结 节 的 结 节 部 分 和 血管 部 分 分 
开 ， 而 区 域 生长 算法 不 能 。 图 4-12 是 分 别 用 本 章 算法 分 割 4 套 粘连 型 血管 
结 节 的 三 维 可 视 化 效果 图 ， 球 状 的 物体 为 结 节 ， 管 状 的 物体 为 血管 。 可 以 
看 出 ， 本 章 算法 能 把 结 节 从 血管 上 分 割 下 来 。 


























a) 粘连 血管 型 结 节 的 第 一 层 影像 b) 粘连 血管 
c) 粘连 血管 型 结 节 的 第 三 层 影像 d) 粘连 血管 











型 结 节 的 第 二 层 影像 
型 结 








本 章 算法 分 割 图 4-12 中 的 4 个 结 节 的 最 佳 带 宽 参 数 选择 范围 与 文献 
[31] 和 第 3 章 中 的 算法 分 割 这 些 结 节 时 带宽 参数 选择 范围 的 对 比 见 
表 4-1。 可 以 看 出 ， 本 章 算法 提出 了 一 个 直接 选择 带宽 参数 的 新 方法 而 不 
同 于 文献 [31] 和 第 3 章 用 迭代 最 优化 的 方法 ， 并 在 带宽 参数 选择 范围 内 
选 出 最 佳 带宽 参数 。 如 对 于 图 4-12a 中 的 粘连 血管 型 结 节 ,文献 [31] 算 
法 的 最 佳 带宽 参数 选择 范围 为 0.5 ~13， 而 第 3 章 的 算法 的 最 佳 带宽 参数 
选择 范围 为 2.75 ~ 9.25， 最 后 用 本 章 算法 则 可 以 直接 求 出 最 佳 带宽 参数 
为 3。 
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9) 


网 4-12 
a), d) 本 章 算法 带宽 为 3 





d) 


b), c) 本 章 算法 带宽 为 3.25 





表 4-1 最 佳 带宽 参数 选择 范围 的 对 比 




















图 号 文献 [31] 算法 第 3 章 的 算法 本 章 算 法 
图 4-12a 0.5~13 2.75 ~9.25 3 

图 4-12b 0.5~13 3.0 ~9.75 3. 25 

Al 4-12¢ 0.5 ~13 3.0 ~ 10. 25 3. 25 

图 4-12d 0.5~13 2.75 ~9.25 3 











BSH 基于 改进 主动 形状 模型 的 
分 割 方法 研究 


5.1 主动 形状 模型 


主动 形状 模型 (ASM, Active Shape Model) 是 由 Coote 等 人 中 提出 
的 。 它 是 建立 在 点 分 布 模型 ( PDM，Point Distributed Model) 的 基础 上 ， 
获取 训练 图 像样 本 的 形状 标记 点 分 布 的 统计 信息 〈 即 平均 形状 ) 和 形状 标 
记 点 的 变化 范围 〈 即 形状 上 的 变形 ) 。 将 平均 形状 覆盖 到 目标 图 像 上 ， 在 
变形 允许 的 范围 内 迭代 并 搜索 目标 边缘 。 在 搜索 的 过 程 中 ， 它 首先 对 一 组 
标 有 特征 点 的 训练 集 图 像 的 形状 和 局 部 灰 度 建 模 ; 然后 不 断 调 节 形 状 和 姿 
态 参 数 从 而 使 最 终 的 形状 达到 最 优化 。 

ASM 包括 形状 模型 、 表 面 模型 和 优化 算法 。 首 先 ， 在 训练 集 样本 图 像 
E, 选取 标记 点 ; 然后 ， 通 过 使 用 旋转 、 缩 放 和 平移 缩小 样本 之 间 的 差 
异 ; 最 后 ， 用 主 成 分 分 析 (PCA, Principal Component Analysis ) ' 引 方法 得 
到 样本 集 的 平均 形状 及 变形 范围 。 以 上 阐述 的 为 建立 形状 模型 步骤 。 表 面 
模型 是 通过 描述 每 一 标记 点 周围 典型 图 像 结构 得 到 的 ， 它 是 以 标记 点 为 中 
心 ， 沿 着 垂直 于 轮廓 的 方向 采样 像素 剖面 得 到 。Cootes 等 人 使 用 这 些 剖 面 
的 归 一 化 的 一 阶 导 数 建立 了 表面 模型 '”; 。 最 后 ，ASM 将 分 割 问 题 转 化 为 
用 优化 算法 如 局 部 搜索 方法 ， 以 及 最 小 均 方差 参数 近似 求 模 型 最 优 解 。 

ASM 的 构成 要 素 有 : 基于 标记 点 的 形状 表示 (采用 标记 点 表示 )、 形 
状 对 准 (采用 广义 Procrustes) 、 基 于 PCA 的 形状 模型 、 由 归 一 化 的 图 
像 梯 度 或 亮度 剖面 生成 的 表面 模型 。 

传统 的 形状 图 形 是 用 手工 标定 在 形状 上 的 标记 点 构成 ， 既 费时 且 容 易 
出 错 。 对 于 三 维 图 像 ， 因 为 需要 分 析 的 数据 量变 大 ， 相 应 地 标记 点 的 数量 
也 惊人 的 大 。 如 文献 [ 36 | 是 用 最 小 描述 长 (MDL，Minimum Description 
Length) 的 方法 建立 标记 点 ， 该 方法 能 自动 地 标记 形状 上 的 标记 点 。 另 一 
种 方法 是 在 三 维 图 像 上 生成 密集 的 三 角形 ， 由 此 选择 标记 点 。 还 可 利用 归 
一 化 互信 息 的 自由 弹性 配 准 技术 在 不 同 的 训练 形状 间 传 播 标 记 点 。 此 外 ， 
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最 小 化 形状 模型 也 可 用 于 选择 标记 点 。Lim 等 建立 了 一 种 三 维 统计 形状 
模型 的 生成 方法 ， 它 首先 从 训练 形状 的 距离 变换 中 建立 平均 形状 ， 接 着 使 
用 四 面体 的 方法 在 平均 形状 中 自动 算出 标记 点 ， 最 后 借助 于 距离 标记 的 方 
法 将 这 些 标记 点 标记 到 每 一 个 训练 形状 上 。 

基于 2- 范式 的 形状 对 准 方法 在 ASM 中 受到 了 广泛 关注 ， 但 该 方法 易 
受到 基于 标记 点 的 影响 。 后 来 有 人 提出 基于 无 穷 - 范 式 或 1- 范 式 的 形状 对 
准 方法 ， 以 便 得 到 更 符合 现实 的 模型 。 

主 成 分 分 析 (PCA) 能 产生 较 好 的 形状 变化 分 解 ， 然 而 这 是 在 假定 误 
差 服从 多 变量 高 斯 分 布 的 前 提 下 。 而 当 假 设 不 满足 时 ，PCA 会 产生 不 合理 
的 形状 分 解 。 

在 ASM 中 ， 通 过 每 一 个 标记 点 并 且 垂 直 于 目标 轮廓 的 纵 训 线段 或 纵 
剖面 上 的 归 一 化 的 图 像 亮 度 或 梯度 均 表示 了 该 点 周围 的 邻 域 特征 ， 而 所 有 
训练 样本 的 对 应 点 的 纵 剖 面 上 的 归 一 化 的 图 像 亮度 或 梯度 经 过 统计 学 习 
后 ， 就 形成 了 该 点 的 局 部 表面 模型 。 

形状 模型 、 局 部 表面 模型 和 搜索 机 制 是 ASM 的 3 个 主要 组 成 部 分 。 
如 果 局 部 表面 模型 是 线性 的 ， 搜 索 方式 是 在 图 像 中 使 用 马 氏 ( Mahalano- 
bis) 距离 测度 将 每 一 个 标记 点 移动 到 其 最 佳 的 匹配 位 置 上 。 当 局 部 表面 
模型 是 非 线 性 的 时 ， 需 要 用 -近邻 分 类 器 以 搜索 标记 点 的 最 佳 位 置 。 通 
常 ,， 为 了 增加 鲁 棒 性 ，ASM 的 搜索 机 制 是 在 多 分 辨 力 框 架 下 实现 。 为 了 使 
ASM 能 够 有 效 地 解释 训练 集 所 蕴含 的 目标 形状 ， 必 须 对 ASM 的 PCA 模型 
进行 恰当 截取 。 


5.2 改进 ASM 在 肺 区 分 割 上 的 应 用 












































5.2.1 肺 区 分 割 研究 现状 


肺 区 分 割 是 肺癌 计算 机 辅助 诊断 的 首要 步骤 ,但 发 生 肿 块 ( 指 直径 大 
于 3cm 的 类 球形 病灶 ) 与 胸壁 粘连 时 ， 会 给 分 割 带 来 困难 。 传 统 方法 大 
多 依据 计算 机 断层 成 像 (CT, Computer Tomography) 值 的 差异 来 分 割 肺 区 
CREEK) 和 胸壁 〈 高 密度 区 域 ) ， 而 这 只 适合 于 分 割 胸壁 无 肿块 粘 
连 的 肺 区 。 如 果 胸 壁 与 肿块 粘连 ， 由 于 两 者 CT 值 接近 ， 会 导致 基于 局 部 
低级 特征 的 传统 分 割 方法 得 不 到 正确 结果 (图 5-1b) 。 近 年 来 ， 很 多 高 级 
的 分 割 方法 应 用 到 医学 图 像 的 分 割 中 ， 如 将 图 理论 、 马 尔 科 夫 随 机 场 应 用 
到 PET-CT 的 肺 肿块 分 割 中 ， 针 对 胸壁 粘连 型 肺 结 节 ( 直径 在 3mm ~ 
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3cm) 的 分 割 方法 已 有 较 多 报道 。 如 文献 [39] 用 形态 学 “ 滚 球 ” 法 填充 
胸壁 粘连 型 肺 结 节 区 域 。 文 献 [40] 利用 肋骨 与 胸壁 之 间 曲 率 一 致 的 特点 
来 分 割 含 胸壁 粘连 型 肺 结 节 的 肺 区 。 文 献 [41] 采用 贝 塞 尔 曲面 将 经 仿 射 
变换 的 平滑 肺 轮廓 模型 拟 合 到 目标 肺 上 ， 再 用 主动 轮廓 算法 加 以 改进 。 但 
该 方法 难以 适用 于 结 节 出 现在 肺 顶 或 肺 底部 的 情形 。 文 献 [42] 先 用 移动 
立方 体 (marching cube) 算法 建立 表面 模型 ， 然 后 用 曲率 识别 并 删除 “有 
问题 ”的 区 域 ， 再 用 径 向 基 洱 数 (RBF, Radiabasis Function) 修复 并 拟 合 
表面 。 但 以 上 方法 只 适用 于 小 的 病灶 区 域 的 填充 ， 难 以 处 理 含 有 大 面积 ; 
KE (如 肿块 ) 的 肺 区 。 文献 [43] 基于 图 集 的 弹性 配 准 来 分 割 含 有 胸壁 
粘连 型 肿块 的 肺 区 ， 分割 准确 率 高 ， 且 能 在 缺少 大 量 正常 肺 组 织 (不 大 于 
30% ) 的 情况 下 维持 完整 的 肺 区 形状 ; 但 处 理 时 间 较 长 ， 在 临床 上 不 能 接 
Zo 文献 [44] 用 肺 影像 库 和 配 准 的 方法 做 含 病 变 的 肺 区 分 割 。 























图 5-1 提取 肺 边缘 标记 点 
a) 不 含 粘 连 胸壁 型 肿块 的 肺 pb) 含有 粘连 胸壁 型 肿块 的 肺 








为 解决 上 述 问题 ， 需 将 先 验 肺 区 形状 模板 导入 到 分 割 模型 中 。 目 前 ， 
基于 先 验 形状 的 分 制 算 法 主要 有 两 类 : 第 一 类 是 带 形状 约束 项 的 水 平 集 
level-set 模型 5]， 第 二 类 是 基于 统计 的 变形 模型 ， 如 主动 形状 模型 
(ASM, Active Shape Model)“|。 前 一 种 方法 因 其 过 于 限定 目标 形状 的 特 
性 而 难以 推广 。ASM 是 一 种 基于 标记 点 分 布 模型 (PDM) 的 形状 表达 模 
型 ， 它 使 用 主 成 分 分 析 (PCA) 学 习 训 练 集 的 平均 形状 和 变形 范围 ， 故 适 
用 于 大 致 轮廓 相似 又 有 一 定 个 体 差异 的 医学 解剖 部 位 分 制 。 但 ASM 分 割 
效果 会 受到 异常 标记 点 (即位 于 伪 边 缘 上 的 标记 点 ) 的 影响 。 异 常 标记 点 
违背 了 ASM 模型 中 元 余 错 误 服从 高 斯 分 布 的 假设 。 因 此 ， 有 效 地 处 理 异 
常 标记 点 是 将 ASM 应 用 于 分 割 肺 区 问题 的 前 提 。 

近年 来 ， 许 多 学 者 在 检测 异常 标记 点 上 做 了 大 量 工 作 。 如 文献 [47] 
用 形状 表达 模型 的 变化 参数 来 检测 标记 点 ， 即 在 形状 变化 参数 下 ， 用 单个 
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标记 点 更 新 的 位 移 量 与 所 有 标记 点 更 新 的 平均 位 移 量 相 比 较 ; 如 果 相 差 很 
大 ， 说 明 该 标记 点 为 异常 标记 点 。 文 献 【48] 用 相 邻 标记 点 间距 离 比 例 作 
为 形状 描述 手段 来 检测 异常 标记 点 ; 如 果 距 离 比例 值 在 训练 出 的 形状 变化 
容忍 模型 的 置信 区 间 外 ， 说 明 该 标记 点 为 异常 标记 点 。 文 献 [49] 将 标记 
点 集 随 机 分 为 若干 个 子 集 ， 若 待 分 制图 像 的 标记 点 子 集 重 建 出 的 形状 与 训练 
样本 集 对 应 标记 点 重建 出 的 形状 差别 较 大 ， 说 明 该 子 集 中 含有 异常 标记 点 。 


5.2.2 提取 标记 点 


首先 ， 用 传统 方法 分 割 肺 区 。 具 体 步 骤 为 迭 代 阅 值 分 割 、 区 域 生长 和 
连通 性 排除 肺 外 区 域 、 剔 除 气 管 、 填 充 内 部 空洞 和 剔除 肝脏 、 左 右 肺 精确 
分 离 分 割 。 

然后 ， 分 别提 取 左 、 石 肺 轮廓 的 标记 点 。 对 于 左 肺 ， 选 择 轮廓 最 上 边 
的 点 、 最 下 边 的 点 、 最 右边 的 点 、 上 半 部 分 最 左边 的 点 和 下 半 部 分 最 左边 
的 点 作为 主 标记 点 。 对 于 右 肺 ， 选 择 轮廓 最 上 边 的 点 、 最 下 边 的 点 、 最 左 
边 的 点 、 上 半 部 分 最 右边 的 点 和 下 半 部 分 最 右边 的 点 作为 主 标记 点 。Aa 
在 两 个 主 标记 点 间 等 距离 间隔 选取 20 个 标记 点 ; 这 样 ， 一 个 轮廓 线 上 及 
=96 个 标记 点 (图 5-1)。 在 图 5-1 中 ， 提 取 肺 边缘 标记 点 中 鸭 出 点 为 主 标 
记 点 ， 亮 点 的 轮廓 线 均 为 其 他 标记 点 。 


5.2.3 建立 形状 模型 


基于 普 鲁 克 分 析 方 法 ( Procrustes Analysis) ， 在 对 章 训练 样本 集 的 各 个 
形状 轮廓 后 ， 求 训练 样本 集 的 平均 形状 过 和 协 方差 矩阵 C; 并 计算 C 的 特征 
值 和 特征 向 量 。 将 特征 值 降序 排列 并 选取 与 前 n 个 最 大 特征 值 分 别 相 互 对 应 
的 特征 向 量 p = [pi ,py，,…,p,] ， 以 及 建立 一 个 线性 可 变 统 计 形 状 模型 x = 元 + 
pbo AP, b 为 形状 参数 。 如 图 5-2 R, HIRR 1 为 训练 样本 集 的 均值 轮 

150 150 150 
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-200 0 200 -200 0 200 -200 0 200 
图 5-2 基于 ASM 模型 求 得 的 变形 形状 
1 一 平均 形状 2 一 变形 形状 
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Rx; 轮廓 线 2 为 基于 该 形状 模型 ， 即 当 分 别 为 1、2、3 时 ， 求 得 的 变形 
形状 。 可 以 看 出 ， 随 着 n 越 大 ， 获 得 的 变形 形状 就 越 和 平均 形状 相似 。 


5.2.4 异常 标记 点 选择 


由 于 肿块 与 胸壁 粘连 ， 基 于 传统 分 割 方 法 获得 的 肺 边缘 轮廓 有 四 口 
(图 5-1b)， 在 四 口 上 的 标记 点 为 异常 标记 点 。 这 些 异常 标记 点 在 ASM 搜 
索 肺 边缘 的 过 程 中 ,会 被 更 新 到 肺 实质 和 肿块 的 过 渡 区 域 。 

首先 ， 将 基于 标记 点 对 齐 含 大 肿块 的 肺 轮廓 Shape, (x, ,yi ，…% ,Yi) 和 
训练 样本 的 均值 轮廓 Shape, (wi ,zi ，,…wi ,zi)。 对 齐 的 步骤 是 ,删除 平移 成 
分 、 缩 放 尺 度 和 旋转 变换 。 

1) 删除 平移 成 分 的 式 为 


(x, 一 > 一 Dyk) (5-1) 
2) 缩放 尺度 的 式 为 
(Cx =x)/s, (y -=y)/s) (5-2) 
式 中 : s =sgt(((z -%)° + (y1 =y)? ++) Ak) 
3) 旋转 变换 中 旋转 角度 的 式 为 
0 = arctan È (win Bi (5-3) 
,1 Wi% + ¥;2;) 


对 齐 的 结果 如 图 5-3 所 示 ， 其 中 标记 点 1 表示 训练 样本 集 的 均值 轮 
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图 5-3 ”对 齐 训练 样本 集 的 均值 轮廓 和 含 大 肿块 的 肺 轮廓 
1 一 均值 轮廓 ”2 一 含 大 肿块 的 肺 轮廓 
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JD, Prick 2 表示 含 大 肿块 的 肺 轮廓 。 

对 齐 之 后 ， 寻 找 在 搜索 过 程 中 可 能 陷 和 人 肿块 区 域 的 异常 标记 点 。 在 含 
大 肿块 的 肺 轮廓 标记 点 集中 寻找 与 均值 轮廓 的 标记 点 Ps (x,y) 距离 最 近 
的 标记 点 ge (w,z) ， 称 为 对 应 标记 点 。 如 果 对 应 标记 点 间 的 距离 大 于 效 
fe (图 5-4) ， 说 明 该 标记 点 为 异常 标记 点 ， 即 


min ,Ix*-gql,>e (5-4) 


x € neighbour( p 

















w 2 3 æ 5 6 70 80 90 100 
标记 点 序列 


图 5-4 训练 样本 集 的 均值 轮廓 与 含 肿块 肺 轮廓 的 对 应 标记 点 距离 





5.2.5 搜索 肺 区 边缘 











首先 ， 需 要 设 定 搜索 范围 即 包 围 核 。 这 是 因为 CT 肺 影像 中 除了 有 肺 
K, DAK, Wass, 包围 核 为 最 上 、 最 下 、 最 左 和 最 右 的 主 标 记 点 
形成 的 矩形 区 域 。 搜 索 的 初始 曲线 为 均值 轮廓 。 当 前 搜索 轮廓 在 包围 核 的 
外 面 时 ， 标 记 点 的 剖面 邻 域内 最 内 部 的 点 为 下 次 迭代 的 位 置 。 

然后 ， 分 段 设 定 搜索 规则 。 对 于 异常 标记 点 ， 搜 索 时 只 需要 考虑 服从 
形变 范围 -3 VA <b< +3 人 ， 且 = 人 (xz -7z)。 其 中 ,图 为 训练 集中 变 
形 的 特征 向 量 和 矩阵; A 为 其 对 应 的 特征 值 ; x 为 训练 集 均匀 轮廓 ; x 为 当 
前 搜索 轮廓 。 对 于 非 异 常 标记 点 ， 除 了 服从 形变 范围 外 ， 还 应 服从 局 部 梯 
度 模型 。 对 标记 点 1 在 剖面 邻 域内 的 每 个 采样 候选 点 ， 计 算 其 梯度 值 L, 
与 标记 点 1 的 局 部 梯度 模型 中 平均 值 1, 之 间 的 马 氏 距离 : 


R= y, Ga) (5-5) 
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并 使 Li) 最 小 的 点 为 下 次 迭代 的 位 置 。 


5.3 数据 集 及 分 割 结 果 


5.3.1 数据 集 


从 广州 医科 大 学 附属 第 一 医院 收集 30 个 含有 粘连 胸壁 型 肿块 的 多 层 
螺旋 CT 胸部 扫描 序列 影像 。 序 列 影像 的 尺寸 为 512 x 512 x (100 ~ 200) 
层 ， 层 厚 为 0.6mm x0.6mm x2mm (或 3mm)。30 个 左 肺 区 中 有 19 个 含 
有 粘连 胸壁 型 肿块 ，11 个 不 含有 粘连 胸壁 型 肿块 。30 个 右 肺 实质 区 域 中 
有 21 个 含有 粘连 胸壁 型 肿块 ，9 个 不 含有 粘连 胸壁 型 肿块 。 左 肺 区 中 肿块 
的 平均 最 大 直径 为 53. 0mm 和 右 肺 区 中 肿块 的 平均 最 大 直径 为 38. 8mm。 
不 含 粘连 胸壁 型 肿块 的 左 、 右 肺 区 ， 作 为 训练 集 样本 分 别 参 与 训练 左 、 右 
肺 轮廓 形状 模型 。 

5.3.2 分 割 结果 

图 5-5 显示 了 基于 本 章 方法 在 肺 实质 边缘 轮廓 搜索 近代 过 程 中 的 初始 
位 置 、 中 间 位 置 和 最 终 位 置 。 图 5-6 对 比 显示 了 改进 ASM 方法 和 传统 方 
法 分 割 带 有 粘连 肺 壁 型 肿块 肺 区 的 效果 。 图 5-7 显示 了 改进 ASM 方法 分 割 
肺 内 含 大 肿块 的 肺 区 的 效果 。 


















































图 $-$ 基于 改进 ASM 的 含 粘连 胸壁 型 肿块 的 肺 区 轮廓 搜索 过 程 
a) 初始 位 置 b) 中 间 位 置 c) 最 终 位 置 
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图 5-6 基于 传统 方法 和 改进 ASM 对 带 有 粘连 肺 壁 型 肿块 的 肺 区 域 分 割 结果 对 比 
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K 5-7 ”基于 改进 ASM 的 对 肺 内 含 大 肿块 的 肺 区 分 割 结 














5.3.3 分 割 结果 评价 与 讨论 
将 分 割 算法 得 到 的 轮廓 包含 区 域 与 医生 手工 绘制 轮廓 包含 区 域 (作为 
“ 金 标 准 ”) 相 比 较 ， 计 算 履 盖 程 度 Overlap, 3 (5-6) 如 下 所 示 : 


Overla 100% 5-6 
p = R „UR, x o (5-6) 


式 中 : Ri 和 RR, 分 别 指 “ 金 标 准 ” 的 像素 集 和 分 害 算 法 的 像素 集 。 
改进 ASM FV AVE REPS EDY 93.6% ， 而 传统 ASM WALI (A 
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为 65. 5% 。 得 到 这 样 的 结果 ， 是 因为 传统 ASM 对 异常 标记 点 敏感 ; 而 本 
章 方 法 能 处 理 异常 标记 点 的 影响 。 

本 章 方法 被 集成 在 基于 ITK 和 VTK 开发 的 医学 影像 后 处 理 的 平台 上 。 
运行 在 配置 为 二 -3210M CPU, 6G 内 存 且 安装 Windows7 操作 系统 的 计算 机 
上 。 实 验 中 ,每 层 平均 分 割 时 间 为 1. 2s， 一 序列 CT 影像 的 分 割 时 间 不 超 


过 2min。 
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